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ABSTRACT 

We analyze several large samples of Active Galactic Nuclei (AGN) in order to establish the 
best tools required to study the evolution of black hole mass (Mbh) and normalized accretion 
rate (L/Lgdd)- The data include spectra from the SDSS, 2QZ and 2SLAQ public surveys at 
z < 2, and a compilation of smaller samples with < z < 5. We critically evaluate the usage of 
the Mg II A 2798 and C IV X 1549 lines, and adjacent continuum bands, as estimators of Mbh 
and i/Ledd' by focusing on sources where one of these lines is observed together with Hj3 . We 
present a new, luminosity-dependent bolometric correction for the monochromatic luminosity 
at 3000A, L3000, which is lower by a factor of ~ 1.75 than those used in previous studies. 
We also re-calibrate the use of L3000 as an indicator for the size of the broad emission line 
region (/?blr) and find that R BLR oc ^3qqq, in agreement with previous results. We find that 
FWHM (Mgll) ~ FWHM (H/3) for all sources with FWHM(Mgll) < 6000 km s _1 . Beyond 
this FWHM, the Mg II line width seems to saturate. The spectral region of the Mg II line can 
thus be used to reproduce H/3-based estimates of Mbh and L/Leaa, with negligible systematic 
differences and a scatter of ~0.3 dex. The width of the C IV line, on the other hand, shows no 
correlation with either that of the H/3 or the Mg II lines and we could not identify the reason 
for this discrepancy. The scatter of Mbh(Civ), relative to Mbh(H/3) is of almost 0.5 dex. 
Moreover, 46% of the sources have FWHM (C iv) < FWHM (H/3 ), in contrast with the basic 
premise of the virial method, which predicts FWHM (C iv) /FWHM (H/3 ) ~ \f3J, based on 
reverberation mapping experiments. This fundamental discrepancy cannot be corrected based 
on the continuum slope or any C iv-related observable. Thus, the C IV line cannot be used to 
obtain precise estimates of Mbh- We conclude by presenting the observed evolution of Mbh 
and L/Ledd with cosmic epoch. The steep rise of L/Lgdd with redshift up to z — 1 flattens 
towards the expected maximal value of L/Lgdd — 1> with lower-MBH sources showing higher 
values of L/L^d at all redshifts. These trends will be further analyzed in a forthcoming paper. 
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1 INTRODUCTION 

The growth of Super-Massive Black Holes (SMBHs), which reside 
in the centers of most large galaxies, progresses through episodes 
of radiatively-efficient accretion, during which such systems ap- 
pear as Active Galactic Nuclei (AGN). Many details of this growth 
process are still unknown. In the local Universe most AGN are 
powered by lower-mass BHs , with Mrh ~ 10 6 ~ 10 8 M,ri, grow- 



NT07). In contrast, the few luminous z ~ 6 QSOs for which Mbh 
was measured show extremely high masses, of about ~ 10 9 M, ; „ 
and high accretion rates, near their Eddington limit. These systems 
could have grown to be the most massive BHs known (Mbh 



ing at very slow rates (e.g., Marconi et al.l 1 20041 : lHasinger et al.l 



l2005l : lNetzer & TrakhtenbrofeOOl hereafter NT07). It is thus clear 
that the more massive BHs accreted at higher rates in the past. 
Indeed, several studies suggest that the typical n ormalized accre- 
tion rates (L/L^dd) increase with redshift (e.g., iFine et al]l2006l 



1Q 10 Mb) as early as z ~ 4 jWillott et al.ll2010l : iKurk et alj|2007t 
|Pe Rosa et al.11201 lMTrakhtenbrot et al.ll201 ll hereafter Tl 1). 

For un-obscured, type-I AGN, Mbh can be estimated by so- 
called "single-epoch" or "virial" estimators. These methods rely on 
the results of reverberation-mapping (RM) experiments, which pro- 
vide empirical relations between the emissivity-weighted radius of 
the Broad Line Region (BLR) and the source luminosity L. These 
^blr ~ L relations are parametrized as ^blr x (^a)"- Assum- 
ing the motion of the BLR gas is virialized, the single-epoch Mbh 
estimators take the general form 
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where Vblr is a probe of the BLR velocity field and / is a fac- 
tor that depends on the geometrical distribution of the BLR gas. 
A common estimator of this type is the "H/3 method" (hereafter 
MbhU^ioo])- Here /?blr is estimated from XL^ at 5100A (here- 
after Ls mo), FWHM(Hg) is the velocity pro xy, and a = 0.6 ± 0. 1 
jKaspi et al.l l200dl . l2005t iBentz et al.l l2009h . Another method is 
based on the Mg II A2798 line and the adjacent continuum luminos- 
ity (L3ooo)- Although the few Mg II-dedicated reverberation cam- 
paigns have not yet revealed a robust R m r — Z/jnoo relation (e.g., 
IClavel et alll99ll ; lMetzroth et alj|2006h . several studies calibrated 
the Mg II relation by using UV spectra of sources that have H/3 re- 
v erberation data. One part icular example is th e relation presented 



McLure & Jarvisl d2002l) . and later refined bv lMcLure & Dunlocl 



d2004 hereafter MP04). The "H/3" and "Mg II" estimators were 
used in numerous papers to derive Mbh for many thousands of 
sources. Th is means focusing on either z < 0.75 or 0.75 < z S, 2 
AGN (e.g., ICroorn et alj|20oi iFine et alj|2006l : IShen et alj|200l ; 
iRafiee & HalftOl ]]). Other studies used these estimators for small 
samples of sources at higher redshifts, where the l ines are ob- 
served in one of the NIR bands dShemmer et al.ll2004l;lNetzer et al.l 
20071: iKurk et alj|2007l ; iMarziani et alj|2009l ; iDietrich et alj|2009t 



Willott et alfcOldT 



In principal, Mbh can also be estimated from the broad 
C IV A 1549 linfl since Z?bl r(C IV) is known from RM experi- 
ments (e^, Kaspi et al .120071). A specifi c calibration of this type is 
given in I Vestergaard & Peterson! J2006I) and several other papers. 
Such methods would potentially enable the study of large sam- 
ples of AGN at 1.5 ^ z 5 for which C IV is observed in large 
optica l surveys (e.g-. IVestergaard et al.l2 008; Vesterga ard & Osmed 
120091) . However, there is evidence that such C IV - based estimates 
are highly uncertain. In particular. iBaskin & LaoJ J2005L hereafter 
BL05) found that the C IV line is often blue-shifted with respect 
to the AGN rest-frame, and that FWHM(ClV) is often smaller 
than FWHM(H/3), both indicating that the dynamics of the ClV- 
em itting gas may be domin a ted by non-virial motion A later study 
by IVestergaard & Peter son (2006) claimed that at least some of 
these findings may be due to the inclusion of narrow-line objects 
and low-quality spectra in the BL05 sample. Despite this reser- 
vation, several subsequent studies of large, flux-limited samples 
clearly demonstrated that the relation between the widths of the 
C IV and Mg II lines shows considerable scatter and, as a result, the 
C IV-based estimates of Mbh can diff er from those deduced from 
Mg II by up to an order of magnitude JShen et alj|2008l :l Fine et al .1 
|2010|) . Moreover, studies o f small samples of z ~ 2 — 3.5 AGN 
by, e.g., jNetzer et al.ll2007l. N07 hereaf ter), JShemmer et"ai]|2004 
S04 hereafter) and ( IDietrich et~al]|2009l . D09 hereafter), show that 
the large discrepancies between C IV and H/3 mass estimates per- 
sist even in high-quality spectra of broad-line AGN (i.e. where 
FWHM(H/3) S i 4000kms~ 1 ). The se issues were investigated in a 
recent study by lAssef et al which suggested an empirical 

correction for the C IV-based estimators that is based on the shape 
of the observed UV-optical spectral energy distribution (SED). This 
correction, however, may turn impractical for large optical surveys 
of high redshift AGN, where only the rest-frame UV regime is ob- 
served. 

In this paper, we re-examine the methods used to derive Mbh 
and L/Ledd of high-redshift type-I AGN. We discuss both large and 
small samples, including some that were never presented in this 
context. The various samples are described in ^2] and the measure- 



Hereafter we refer to the UV lines under study simply as Mg II and C IV. 



ment procedures in 5j3] We discuss the monochromatic luminosities 
and bolometric corrections in ij4] We briefly discuss the premise of 
the virial assumption in for estimating Mbh in m SHwe exam- 
ine how the Mg II emission line complex can be used to measure 
Mbh, an d in S we provide new evidence regarding the fundamen- 
tal limitations of the C IV method. Finally, in ij8]we briefly describe 
the observed evolution of Mbh and L/L^dd, as measured with the 
tools developed in this paper, and summarize our conclusions. A 
more detailed analysis of the evolutionary trends is deferred to a 
forthcoming paper. Throughout this work we assume a cosmologi- 
cal model with fl A = 0.7, Cl M = 0.3, and H = 70km s -1 Mpc~'. 



2 SAMPLES SELECTION 

The main goal of the present work is to test how the Mgll A2798 
and Civ A 1549 emission line complexes can be used to estimate 
Mbh, ^bol and L/L^dd, and to apply these methods to probe the 
evolution of L/L^dd to z — 2. While Mg II can be compared to H/3 
or Civ within the same spectroscopic window (at z ~ 0.6 or z ~ 
1.7, respectively), the comparison between C IV and H/3 has to be 
based on a combination of separate observations. This dictates two 
distinct types of samples: (1 ) large samples d rawn from the Sloan 
Digital Sky Survey (SPSS: (York et al.ll200d) . 2dF QSO Redshift 
survey (20Z: ICroom et alj2004h and t he 2dF SPSS LRG And QSO 
survey (2SLAO: lRichards et al]|2005l : ICroom et alj200"8T) catalogs; 



and (2) several smaller AGN samples, with publicly available data. 
The following details all the samples used in the paper, which are 
summarized in TableQ] We note that the 2QZ and 2SLAQ surveys 
provided relatively few high-quality spectra, which are usable for 
the comparisons between the different lines (see below). All the 
spectra analyzed in the pres ent work were correc ted for galactic 
ex tinction using the m aps of ISchlegel et al .1 i ll 9981) and the model 
of lCardellietall ( ll989l) . 



2.1 SDSS DR7 

We queried t he public catalogs of t he seventh data release of the 
SPSS (PR7; lAbazaiianetal.ll2009h for all sources that are clas- 
sified as "QSO" and have a high confidence redshift determina- 
tion (i.e. zconf > 0.7) in the range of z < 2. This resulted in 
74517 sources. First, we omitted from this sample 7845 objects 
for which the line fitting procedures resulted in poor-quality fits, 
similarly to our criteria in NT07 (see more details in 5j3}. Next, 
we omitted from our analysis all radio-loud sources. The radio 



loud ness of each source, R, = f v (5 GHz) / f v (4400A) (follow- 
ing Kenennam^TalJ ^8^_was derived from the FIRST data 



jBecker et alJll995l : IWhite et al.lll997l) incorporated in the SPSS 
public archive. For sources at z > 1, we estimated f v (4400A) 
using the continuum flux density near 2200A and by assuming 
f v oc v -°- 44 dVanden Berk et all 1200 ll VdBOl hereafter). We re- 
moved all sources with i?L > 10 from our sample. 

Finally, we omitted 1764 sources that are classified as broad 
absorption line QSOs (BALQSOs), based on their Mgll and/or 
Civ spectral regions. Fo r this, we used t he relevant flags in 
the SPSS/PR7 catalog of IShen et all ( feOlll. see below) that, in 
turn, is large l y base d on the SPSS/PR5 catalog of BALQSOs of 
iGibson et alj J2009h . with further (manual and hence partia l ) clas- 
sification of the post-PR5 spectra. Although the lShen et alj fcoill) 
BALQSO classification is probably incomplete, its combination 
with our fit quality criteria provide a sample which is almost com- 
pletely BALQSO-free. 
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The final SDSS sample comprises 20894 for which Hj3 is ob- 
served (z < 0.75) and 43995 sources for which only the Mg II line 
is observed (i.e. 0.75 < z < 2). There are 6731 sources for which 
both the H/3 and Mg II lines are observed (0.5 < z < 0.75; the SDSS 
"Hj3xMgn" sub-sample hereafter) and 10910 sources where both 
the Mgll and ClV lines are observed (1.5 < z < 1.95; the SDSS 
"Mg II xC IV" sub-sample hereafter). 

The gene ral SDSS sample s ignificantly overlaps with the sam- 
ple studied bv lShen et aildioT lh . Virtu ally all our z > 0.75 sources 
are found within the lShen et~ U2011I) catalog, in particular about 
98% of our SDSS "H/3xMgH" and "MgllxClV" sub-samples, 
which have Mbh derived from more than one emission line in that 
catalog. At lower redshifts the fractions are lower, and o nly ~ 68% 
of our Hj8-only SDSS sources can be found within the lShen et al .1 
d201 lh catalog. This is due to our choi ce to includ e relati vely faint 
sources that were excluded from the IShenetal.1 (1201 lh catalog , 
namely sources with M; > —22 (see also Schneider et al.l [2oTof> . 
We verified that these faint sources have high-quality Hj8 fits, ac- 
cording to our criteria (see i}3j. In addition, we note that our SDSS 
"H/3 xMg II" sub-sampl e includes 354 of t he 495 sources (~ 75%) 
analyzed in the study by Wang et al. I J20091) . The discrepancy is due 
to the fact that lWang et al .1 J2009h chose to include in their sample 
sources with slightly lower redshifts (0.45 < z < 0.75). Our general 
(Hj8) SDSS sample includes 473 of their sources (~ 95%). 



2.2 2QZ 

We used the public 2QZ catalog of ICroom et alj J2004I) . initially 
choosing all 0.5 < z < 2 sources which are classified as "QSO" and 
have high confidence redshift determination (i.e. ZQUAL flag better 
than 22). Our query naturally omits BALQSOs, which are classified 
as such within the 2QZ catalog (b ased on visual insp ection of the 
spectra, which is incomplete; see ICroom et al. I l2004h . This query 
provided 20315 sources. After applying the fit quality criteria to 
this large Mg II 2QZ sample, it comprises 8873 sources. Gener- 
ally, spectra obtained as part of the 2QZ and 2SLAQ AAO-based 
surveys are not flux calibrated. Previous studies used single-band 
magnitudes, in either the Bj or the g bands, to derive monochro- 
matic luminosities, assuming al l spectra follow a uniform UV-to- 
optical SEP of fy ~ y -°- 5 (e.g.. lCroom et alj2004lRichards et all 
l2005l : lFine et~al]|2008h . Here we employed a more robust flux cal- 
ibration scheme, which relies on all the photometric data available 
in the 2QZ and 2SLAQ catalogs. This procedure is discussed in 
Appendix [A] 

In order to expand the luminosity range of the "Hj8 xMg II" 
sample, we preformed a separate query of all 0.35 < z < 0.55 
2QZ sources, ignoring their class and ZQUAL. For the selection 
of "MgllxClV" sources we focused on sources at 1.5 < z < 1.8 
(within the larger Mg II sample) These choices were made in or- 
der to include sources which might be flagged as problematic due 
to the proximity of the Hj8 line to the telluric features near 7000A 
or to the limit of the observed spectral range, and to avoid these 
problems from affecting the Civ line. The "H/3 x Mgll" sources 
which passed the fit quality criteria were manually inspected, to 
omit any non-QSO objects or otherwise unreliable spectra. These 
procedures resulted in merely 238 sources with reliable H/3 fits and 
489 sources with reliable C IV fits. 

To summarize, we have 9111 2QZ sources in our final sam- 
ple, 238 of which have both H/3 and Mg II, and 489 sources with 
both Mg II and C IV. These sources span the magnitude range 
bj ~ 18 — 20.85, which corresponds to the luminosity range 44 < 



log (L bol /ergs _1 ) < 46 atz~ 1.5 (see fl4~2t . The 2QZ sample thus 
reaches about a factor 3.75 deeper than the SDSS sample. 

2.3 2SLAQ 

Sources from the 2SLAQ survey were selected by u sing the pub- 
lic 2S LAQ QSO catalog, compi led and presented in lCroom et al .1 
(2008), and further analyzed in iFine et al. Our query fo- 

cused on 0.5 < z < 1.8 sources which are classified as "QSO", and 
resulted in 5873 sources. Here too, BALQSOs are largely omitted, 
due to their separate classification. Low-redshift 2SLAQ sources 
suffer from the same limitations affecting the 2QZ low-redshift 
sources. We preformed a separate search to locate "H/3 x Mgll" 
candidates at 0.4 < z < 0.6. Our careful manual inspection re- 
sulted in only 39 sources with reasonable fits to both the H/3 and 
Mgll lines. We chose not to include any of these sources in the 
present analysis, due to their poor S/N and small number. Sources 
with 1.5 < z < 1.65 in the large Mgll sample have both Mgll 
and C IV. We have 217 such sources, out of a total of 2287. These 
sources span the magnitude range g ~ 18 — 22, which corresponds 
to 44 < log (Lb i/ergs _1 ) < 46 at z — 1.5. This is almost an order 
of magnitude deeper than SDSS, and about a factor of 2.5 deeper 
than 2QZ. 



2.4 Additional Samples 

There are many smaller samples in the literature with measure- 
ments of more than one of the three emission lines discussed here. 
We choose to focus on samples for which the relevant line mea- 
surements are either preformed by procedures very similar to ours 
and publicly tabulated, or on samples for which such measure- 
ments could be preformed on publicly available spectra. In partic- 
ular, we used the samples of BL05 (59 sources), N07 (and S04; 44 
sources), Shang et al] (2007, hereafter Sh07 hereafter ; 22 sources), 
ISulentic et al.l J2007L hereafter Sul07 ; 90 sources) iMcGill et al.l 
d2008l , hereafter Mc08; 19 sources). Ipietrich et all 12009I , here- 
after D09; 10 sources), iMarziani et alj J20091 hereafter M09; 30 
sources), and Til (40 sources). Basic information on these sam- 
ples is given in Table [T] while Appendix |B1 provides more details 
regarding the acquisition and contents of the samples. In particular, 
Appendix iBldetails which line measurements we used cases where 
sources appeared in more than one sample. 

These small samples probe a wide range of redshift and lumi- 
nosity and most of them have higher S/N than those typically ob- 
tained within the large SDSS, 2QZ and 2SLAQ surveys. In many 
of these samples the different lines were not observed simultane- 
ously. Thus, it is possible that line and/or continuum variability 
contributes to the scatter in the relationships discussed in this work. 



3 LINE AND CONTINUUM FITTING 

We have developed a set of line fitting procedures to derive the 
physical parameters related to the H/3, Mg II and C IV lines. These 
procedures are similar in many ways to those described in previ- 
ous publications (e.g., S04, N07, NT07, S08, F08). In particular, 
the H/3 and Mg II fitting procedures follow those used in NT07 and 
Til, respectively. In short, a linear continuum model is fitted to the 
flux densities in two 20A-wide "bands" on both side of the rele- 
vant emission line. The bands are centered around 4435 & 5110A 
for Hj8, 2655 & 3020A for Mg II and 1450 & 1690A for C IV. For 
the H/3 and Mg II lines, we subtract a template of Fe II and Fe III 
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Table 1. Summary of Samples 



Name (Acc.) 


Nlbtal 


z 


Hj8 


Lines used a 
Mgll 


Civ 


Sizes of cross-calibration samples 
Hj3xMgII MgllxClV Hj3xClV 


SDSS 


66653 


0- 


2 s 


s 


s 


6731 


10910 


2QZ 


9111 


0.35- 


-1.7 s 


s 


s 


238 


489 


2SLAQ 


2287 


0.45- 


1.6 s 


s 


s 


39 


217 



Baskin & Laor (2005. BL05) 


59 


0.25-0.5 


s 




t 






81 


Netzer et al. (2007. N07) b 


44 


2.2-3.5 


s 




s 






44 


Shang et al. (2007. Sh07) 


22 


0.1-0.4 


t 


t 


t 


22 


22 


22 


Sulentic et al. (2007, Sul07) 


90 


< 1 


t/s 


s 


t 


9 


1 1 


57 


McGill et al. (2008, Mc08) 


19 


~0.36 


t 


1 




19 






Dietrich et al. (2009. D09) 


10 


1-2.2 


t 


1 


t 


7 


7 


9 


Marziani et al. (2009, M09) 


30 


1-2.4 


t 


S 


s 


16 


12 


6 


Trakhtenbrot et al. (2011, Till 


40 


~4.8 




s 


s 




40 





a Type of data used: "t" - tabulated data; "s" - spectra fitted by our team (i.e., this work,Jjhemmej^rtal ii 2004, N07 and Til). 
b The majority of sources in this sample were presented in Shemmer et al. (2004). 



emission lines. This is done by choosing the best-fit broadened 
and sh ifted template based on data published bv lBoroson & Greet] 
( 1 19921) and lVestergaardl d2004h . The grid spans a broad range both 
in the width and (velocity) shift of the iron features. In the case 
of Mg II, we add flux in the wavelength range under the emission 
line itself (see Appendix IC2b The main emission lines are mod- 
eled by a combination of two broad and one narrow Gaussians. In 
the cases of Mg II and C IV this model is duplicated for each of 
the doublet components. In Appendix [C] we discuss several minor 
but important technical issues concerning the fitting of the lines, 
including: (1) the ranges of all relevant parameters; (2) the nar- 
row ( lOOOkms -1 ) components of all lines; (3) the improved 
iron template for the Mg II spectral complex and its scaling; (4) 
the difference between our fitted L3000 an d the "real" power-law 
continuum, and (5) the emission features surrounding the C IV line. 

All fitting procedures were tested on several hundred spec- 
tra, selected from all samples used here. For these, we individually 
inspected the fitted spectra. In particular, we verified that the sim- 
plifications involved in the C IV fitting procedure (Appendix IC3t 
are justified. After fitting the SDSS, 2QZ and 2SLAQ samples, we 
filtered out low-quality fits by imposing the criteria ■/} < 5 and 
R 2 > 0.20 Both these quantities were calculated over a range of 
14000kms _1 centered on the peaks of the relevant emission lines, 
and adjusted for degrees of freedom. This choice is similar to the 
one we made in our previous analysis of large samples (NT07). 
The smaller samples were treated manually, securing high-quality 
fitting results for all sources where adequate spectral coverage of 
the various lines was available. We stress that even when two of 
the three main lines are observed in the same spectrum (i.e. for 
the "Hj3xMgH" and "MgllxClV" sub-samples), the fitting pro- 
cedures are ran separately for each emission line complex. Thus, 
the parameters of the different models are fully independent, even 
when it is possible to assume otherwise. 

The line fitting procedures provide monochromatic luminosi- 
ties (L1450, £3000 and/or L5100) an d emission lines widths. For the 
latter, we calculated the FWHM, the line dispersion o~ (following 
IPeterson et all2004T) and the inter-percentile velocity (IPV; follow- 



2 R 2 is the coefficient of determination, related to the amount of variance 
in the data that is accounted for by the model. 



ing lCroom et alj|2004l) . We estimated Relr for all the sources for 
which the "H/3 method" is applicable, by: 



^blr(Hj8)= 538.2 



-5100 



It — days . 



(2) 



10 46 ergs- 

Thisis an updated version of the correlation found by [Kaspi et all 
d2005h. taking i nto ac count the improved measurements presented 
in iBentz et al. I d2009l) , but fitting the RM measurements of only 
the sources with L5100 ^ lO^ergs -1 , typical of the luminosities 
of high-redshift sources. Mrh(H]3) is derived by assuming / = 1 
dNetzer & Marziariill201Ct and references therein) and is given by: 



Mbh (Hj3 ) = 1.05 x 10 s 



-5100 



10 46 



ergs 



0.65 



FWHM(H J 6) 



10 3 kms~ 



M . 
(3) 

We have also experimented with flatter ^blr-^5100 relations, of 
j?BLR ^sioo' as su gfi es ted by some stu dies (e.g. IBentz et all 
|2009|) . We found, as in iKaspi et alj d2005l) , that these have little 
effect on our main results, provided that the constant in Eq.[2]is 
adjusted accordingly. In particular, for the median luminosity of 
the SDSS "HjS xMg II " sub-sample (L^nn = 6 x 10 44 ergs" 1 ) the 
^BLR-isioo relation of IBentz et alj d2009h would have resulted in 
^blr estimates that are larger by 0.03 dex than the ones obtained 
by Eq.f2] For the "HjSxClV" and "Mgllx ClV" sub-samples we 



assumed the reverberation-based relation of Ka spi et al .1 d200' 



tples 
007h : 



^blr(Civ) = 107.2 



It — days . 



(4) 



10 46 ergs- 

We al so experimented with the relation of Ivestergaar d & Peterson! 



BLR x ^1450 > an d verified that none of the results that fol- 



d2006|) . R 
low depend on this choice 



4 LUMINOSITIES AND BOLOMETRIC CORRECTIONS 
4.1 Monochromatic Luminosities and SED Shapes 



The UV-optical continuum emission of type-I AGN, between 
~1300 and 5000A is well described by f v ^ V -0 ' 5 (e.g., VdBOl). 
The typical luminosity scaling values in our SDSS sample are 
{£3000/^5100) = 1-62 and {Z. 145 o/L3ooo) = 1-39, with standard de- 
viations of 0.12 and 0.14 dex, respectively. Assuming the average 
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48^ 



-I 




44- B 3a^''y 



43 ;V 

43 44 45 46 47 48 

^ioo 76 ^ 8 " 1 

Figure 1. A comparison between Z4450 and L5100 for for the different 
"Hj3xClV" samples: BL05 (black squares), N07+S04 (magenta squares), 
Sh07 (red circles), Sul()7 (blue crosses) D09 (magenta crosses) and M09 
(magenta triangles; see Table [T] for references and details). The dashed 
line represents the 1 : 1 relation, and the dot-dashed line represents Z4450 = 
1.97Z.5ioo, the median value found for the N07+S04 sample, and typical 
for o ther high-luminosity sour ces. A power-law continuum with a v = —0.5 
(e.g. lVanden Berk et al.l200lh would imply L1450 = 1.88Z,5ioo- 

conversion of L3000 to the "real" continuum (see Appendix |C2| >, 
these ratios become (L3000/L5100) = 1-41 and (^1450/^3000) = 
1.59, respectively. For comparison, the ratios implied by an a v = 
—0.5 power law are 1.30 and 1.44, respectively. Figure[T]shows that 
the smaller "H/3 x C IV" sub-samples generally follow the scalings 
of a v — —0.5. However, the ratio of ^1450/^5100 — 1-95 provides a 
somewhat better scaling for these sources, for which both L5100 and 
Z4450 are directly observed. This ratio is consistent with the VdBOl 
result, which is based on a composite of sources for which either 
^5100 or ^-1450 were observed. In particular, all the PG quasars tha t 
are part of the reverberation-mapped sample o f iKaspi et al.l fcOOCh 
follow these luminosity scalings (these are part of the BL05, Sh07 
and Sul07 samples). Thus, we do not expect that SED differences 
would play a major role in scaling any relation that is consistent 
with the RM experiments. This is a crucial point in single-epoch 
Mgn determinations and is further discussed in 5j7] 

We note that 87% of the sources in the SDSS "MgllxClV" 
sub-sample have L1450 > L3000, and 93% of the sources in the 
smaller "MgllxClV" sub-samples (Sh07, BL05, N07+S04, M09 
& D09) have I4450 > L 51 oo (i.e., Ct v uyp > —1). In contrast, half 
of the sources (6/12) presented bv lAssef et al.l J201ll) have Z4450 < 
.£,5100- Such extreme UV-optical SEDs may represent AGN with in- 
trinsically different properties or, perhaps, are affected by a higher- 
than-typical reddening. These issues are further discussed in ^7] 

4.2 Bolometric Corrections 

To determine the bolometric luminosity (Lbd) one has to as- 
sign a bolometric correction factor, /bd W = LbolM-^A- Here 
we focus on fi, \ (3000A). Sever al earlier studies assumed a 
constant f hn] f3000A] , e.g., ~ 5.9 felvis et al.l [l994h . or ~ 5.15 
dRichards et al.l2006al) . There are two problems with this approach: 
(1) The constant / bol (3000A) was based on the total, observed X- 



ray to mid-IR (MIR) SED of AGN. As such, it includes double- 
counting of part of the AGN radiation (the MIR flux originates 
from re-processing of the UV-optical radi ation). This results i n 
overestimation of Lb i and thus fi, Q \ (e.g., iMarconi et al.l |2004|) . 
(2) The shape of the SED is known to be luminosity-dependent. 
This dependence is most significan t at the X-ray regim e, and per- 
haps also at UV wavelengths (e.g., IVignali et al.ll2003l) . The gen- 
eral trend is of a decreasing Lx with increasing Luy_ ODt . Thus, 
/~hnl sh ould decrease with increasing UV luminosity. IMarconi et al.l 
(2004) provides a luminosity-dependent prescription for estimating 
/bol(4400A) that addresses these two issues. The prescription can 
be modified to other wavelengths, by assuming a certain UV-optical 
SED. 

We derived a new prescription for fi, \ (3000 A) by using the 
SDSS "Hj8 xMg II" sub-sample. First, we converte d the measured 
L5100 of each source to L(,ol using the prescription of lMarconi et al.l 
d2004l) and assuming f v oc v . This provides, for each of the 
6731 sources, L3000 and ^bol^3000. which are presented in Fig- 
ure |2] The derived bolometric corrections are systematically lower 
than the aforementioned fixed values. For example, the typical 
correction is just 3.4 for sources with log (Z^ooo/ergs -1 ) — 3 x 
10 45 ergs -1 , which is the median luminosity of the SDSS sources 
at z— 1.2. This is a factor of 1.5 lower than the val ue of 5.15 used i n 
several other studies of L/Lgdd at z ~ 1 — 2 (e.g. jFine et alj|20"08r) . 
Fig. [2] further reveals that, despite the large scatter at low L3000, 
there is a clear trend of decreasing ft, a i (3000A) with increasing 
L3000, as expected. To quantify this trend, we grouped the data in 
bins of 0.2 dex in L3000 and assumed that the error on /b i (3000A) 
is (j;/ \/N[, for the i-th bin in Ljooo- An ordinary least squares (OLS) 
fit to the binned data points gives: 

/ b0] (3000A) =4.12- 2. 13^,000,45 (5) 

+ 1.76-^3() 00 ,4 5 _0 - 54 ^3ooo,45 ' 

where ^f 30 oo.45 = log (L 30 oo/10 45 ergs" 1 ) . As Fig. [2] clearly 
shows, this relation predicts relatively large bolometric corrections 
for low-luminosity sources, to the extreme of /b i (3000A) ^ 9.5 
for L3000 % 10 44 erg s~ 1 . Similarly high val ues were reported in th e 
past, for a small minority of sources (e.g., iRichards et alj l2006al) . 
We suspect that the high values predicted by Eg . [5] are the result 
of unrealistic extrapolation of the IMarconi et al.l d2004h Lb — L^ \ 
relation towards low luminosities, where the number of measured 
points is small. Moreover, we note that low-luminosity sources may 
suffer from more significant host-galaxy contamination, which re- 
sults in a systematic over-estimation of L5100 and Lb i- We thus 
caution that the bolometric corrections we provide here should not 
be used for sources with L3000 10 44 ergs _1 , where the emission 
from the host galaxy is comparable to that of the AGN and in partic- 
ular in cases where the spectroscopic aperture affect s the determi- 
nation of the host emission (see, e.g., the analysis in lStern & LaoJ 
120121) . 

Figure [3] compares the bolometric luminosities derived from 
L3000 by using Eq.|5]to those calculated from Lsioo- The two meth- 
ods provide consistent estimates of Lb i with a scatter of less than 
0.09 dex (standard deviation of residuals). The small scatter is prob- 
ably dominated by the range of UV-to-optical slopes of individual 
sources. 

We finally note that the real uncertainties on such estimates of 
Lbol are actually governed by the range of global SED variations 
between sources, as well as the assumed physical (or empirical) 
model for the UV SED. For example, the assumed exponent of the 
X-ray model and the Lx — Ljjy relation may amount up to ~0.2 dex 
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Figure 2. The relation between the empirically derived bolometric cor- 
rections (/ bo i (3000A) = Lbol/i300o) and L3000, for the SDSS "H/J xMg II" 
sub-sample. Large symbols and error bars represent the binned data and the 
corresponding scatter (standard deviations). The solid black line represents 
the best-fit cubic polynomial to the binned data, given in Eq.|5] The dashed 
horizontal line repr esents fb i (3000A) = 5.15, the value used in several 
earlier studies (e.g.. Richards et al. 2006a). 



in Lx, and thus in t he calculated L(, l (see, e.g., IVienali et alj2003l : 
iBianchi et ai]|2009l) . Two very rece nt studies furth e r dem onstrated 
these complications. The study by | Runnoe et al .1 J2012I) showed 
that even the uniform lElvis et all Jl994h and lRichards et al] d2006al) 
SEDs may provide fi, \ (3000A) as low as ~ 3, given that the 
integration is limited to X < 1 flm (thu s neglecting re-proc essed 
emission. However, the best-fit trends of iRunnoe et pre- 
dict /boi (3000A) ~ 5.3 for sources with L3000 ~ 10 4S , c onsistent 
with the commonly used value of 5.15. |jin et al. used an 

accretion disk fitting method that results in a much higher (un- 
observed) far-UV luminosity for a given optical and/or near-UV 
luminosity. Naturally, this model produced very high bolometric 
corrections, with /t, i (5100A) ^ 12 (and as high as ~ 20-30) for 
several sources with L5100 — 3 x 10 44 ergs~', compared to the 
iMarconi et ail d2004h prediction of / bol (5100A) ~ 7. This com- 
pletely new approach to the estimate of Lt, i in AGN will not be 
further discussed in the present paper. Instead, we advice the us- 
age of the cor rections given by Eq. 5] which supplement the optical 
corrections of lMarconi et al] d2004 ). 



5 VIRIAL M BH ESTIMATES: BASIC CONSIDERATIONS 

A main goal of this paper is to present a critical evaluation of 
the various ways to measure Mbh by using the RM-based "virial" 
method. It is therefore important to review the basic premise of the 
method and the justifications for its use. 

Four critical points should be considered: 

(i) The emissivity weighted radius of the BLR, ^blr> ls known 
from direct RM-measurements almost exclusively for only the H/3 
and C IV emission lines, and to a much lesser extent for Mg II (see 
JTJ. The expressions chosen for the present work are given in Eqs.|2] 
and|4] They depend on the measured L5100 and Z4450, and involve 



Figure 3. A comparison of estimates of Lbol based on L3000 and L5100, for 
the SDSS "H/3 xMg II" sample. Large symbols and error bars represent the 
binned data and the corresponding scatter (standard deviations). The dashed 
line represents the the 1 : 1 relation. 



the assumption that the derived luminosities require no reddening- 
related, or other, corrections. The slope a of these correlations 
(^blr x L a ) is empirically determined to be in the range 0.5-0.7, 
by a simple regression analysis of the observational results. It is not 
known whether the fundamental dependence is on LsioOi ^1450' me 
ionizing luminosity or perhaps L(,ol> although there are theoretical 
justifications for all these cases. It is also not clear whether the slope 
itself depends on luminosity (i.e. whether the same relation holds 
for al l luminosities, see, e.g. jBentz et a l. 2009; Netze r~& Marzianil 
|2010|) . Thus, the approach adopted here is to use all quantities as 
measured. 

(ii) The mean ratio of the RM-based H/3 and C IV radii, assum- 
ing the Rblr ~ L relations, is about 3.7. The number depends on the 
measured lags and the mean luminosities of the objects in the RM 
samples used to derive the relations, and may signifi cantly vary for 
indiv idual sources (see, e.g., the range of ratios in lPeterson et al] 
12004b . The mean ratio of the two involved luminosities in the two 
RM samples is L1450/L5100 ~ 1-9. This is similar to the mean ratio 
in several much larger samples where the two wavelength regions 
are observed (e.g., S04, BL05, N07; see 94. It and is also similar to 
the SDSS-based composite spectrum of VdBOl. The working as- 
sumption, therefore, is that this ratio represents the population of 
un-reddened type-I AGN well. Given this, the virial method cannot 
be applied directly to sources where ^1450/^5 100 deviates signifi- 
cantly from this typical value, since in such sources ^blr ma y not 
scale with the luminosity in the same way as in the two RM sam- 
ples used to derive the equations. For example, L1450/L5 100 << 1-9 
may indicate significant continuum reddening, which will result in 
a systematic underestimation of Rblr, if one uses Eq.|4]or similar 
relations. Correcting for such reddening must be performed prior 
to the estimation of /?blr or ^bh- 

(iii) The virial motion of the BLR, combined with the adop- 
tion of line FWHMs as the gas velocity indicator, and the assump- 
tion that the global geometry of the H/3 and C IV parts of the 
BLR are similar (e.g. two spheres of different radii) lead to the 
prediction that FWHM(Civ) /FWHM (H/3) ~ VJ7J. This is con- 
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firmed in a small number of intermed iate luminosity objec ts show- 
ing the expected ^1450/^5100 (e.g.. IPeterson et al]|2004) . Given 
this, the simplest way to proceed to measure Mbh in large samples, 
which lack any additional information regarding the ^1450/^5100 or 
FWHM(Civ) /FWHM(HjS) ratios, is to assume the same is true 
for all sources. 

(iv) The best value of the geometrical factor / in the mass equa- 
tion (Eq. [TJ is 1.0. This is an average value obtained from the 
comparison of RM-based "virial products" and the Mbh ~ O* rela- 
tion, for about 30 low redshift t ype-I AGN (e.g JOnken et alj2004l 
IWoo et alj|2O10l : lGraham eTalfeOl ll) . None of these sources show 
a large deviation from this value. These calibrations rely predomi- 
nantly on H/3, while Civ-related observables contribute to the es- 
timation of virial products in only 5 sources. Unfortunately, there 
are no direct estimates of / that are based solely on C IV. 

Given points (ii) and (iii) above, the value of / used in Mbh 
estimates must be the same for the two emission lines. This is not 
meant to imply that certain sources cannot have two different virial- 
ized regions, for Hj8 and C IV, with different geometries and values 
of /. It only means that the method is based on certain samples 
with certain properties and hence should not be applied to objects 
with different properties. In objects where / is not directly mea- 
sured and the line width ratio deviate significantly from the above, 
e.g. objects with FWHM(H/3) > FWHM(Civ), at least one of the 
lines should not be used as a Mbh indicator within the framework 
of the virial method. Since estimates of / based on FWHM(H/3) 
are known to be correct in most of the H/3-RM sample, we prefer 
to adopt in such cases the assumption that FWHM(H/3) provides 
the more reliable mass indicator. 

The above considerations suggest that a single epoch mass 
determination based on H/3 and L5100 is a reliable mass estimate 
provided there is no significant continuum reddening. In case the 
amount of reddening is known, it should be taken into account prior 
to the application of the method. In the absence of direct mass cal- 
ibration based on lines other than H/3, there are only two alter- 
natives: use theoretical conjectures, or a-priori knowledge about 
the line width. The first can be used for the Mgll line that is 
thought to originate from the same part of the BLR as H/3. If 
FWHM(Mgll) ~ FWHM(H/3), then the line can be used for es- 
timating Mbh- The second can be used for the Civ line since 
^blr(Civ) /-Kblr(HjS) is known. As explained, the requirement is 
FWHM(H J 6) /FWHM(Civ) ~ blr (C iv) /R B lr (H/3) . If this 
ratio is indeed found in the majority of objects with measured 
FWHM(H/3) and FWHM(C IV), then a single epoch mass estimate 
based on C IV can be safely used. In the following sections we use 
such considerations to assess the validity of the use of the Mg II and 
C IV lines as mass indicators in type-I AGN. 



6 ESTIMATING M BH WITH Mg II A.2798 

As discussed in ij5] the fact that there is no RM-based determi- 
nation of i?BLR(Mg II) means that the only way to obtain a Mg II- 
based estimator for Mbh ls to calibrate it against Mbh(H/3), based 
on the the assumption that ^BLR(Mgll) = ^blr(H/3) (verified by 
photoionization calculations). We therefore have to show that L3000 
can be used to estimate ^blr(H/3), and that FWHM(Mgll) ~ 
FWHM(H/3). In what follows, we discuss these relations sepa- 
rately, and evaluate the ability of the (combined) virial product to 
reproduce Mbh(H/3). 



6.1 Using L3000 to Estimate Rblr 

In Figure [4] we present the relation between the calculated val- 
ues of ^blr(H/3) and L3000, f° r all the sources in the different 
"HjSxMgll" sub-samples. There is a clear and highly significant 
correlation between these two quantities. Since Rblr ls calculated 
directly from £5100- this relation reflects the narrow distribution of 
UV-optical continuum slopes. To quantify the ^blr — ^3000 rela- 
tion, we bin the SDSS and 2QZ "HjSxMgll" sub-samples (sepa- 
rately) in bins of 0.2 dex in Z/jooo- The uncertainties on each binned 
data-point are assumed to be the standard deviations of the values 
included in the respective bin. The typical uncertainty is of 0.13 
dex . This is a conservative choice, which attempts to account for 
the entire scatter in the dataQ All other "H/3 xMgll" sub-samples 
remain un-binned. For these, we assume uncertainties of 0. 1 dex on 
^blr and 0.05 dex on £3000- These choices reflect the absolute flux 
calibration uncertainties and the uncertainties related to the contin- 
uum and line fitting processes. Sin ce the uncertainties in both axes 
are com parable, we u se the BCES jAkritas & Ber shadv 199^) and 
FITEXY JPressll2002T) fitting methods, both designed to also take 
into account the scatter in the data. All the BCES correlations are 
tested by a bootstrapping procedure with 1000 realizations of the 
data under study. We used the more sophisticated ver sion of the 
FITEXY method, as presented bv lTremaine et alj d2002l) . where the 
error uncertainties on the data are scaled in order to account for the 
scatter. Thus, all our FITEXY correlations resulted in Xv — 1- The 
best-fit linear relations (also shown in Fig.f4). parametrized as 

log^ BLR = a log ( 10 44 3 e r g S -i ) + P ' (6) 

resulted in a ~ 0.615 and /3 ~ 1.33 for both the BCES and FITEXY 
methods, The exact values and associated uncertainties are given 
in Table [2] The standard deviation of the residuals is about 0. 1 
dex. We note that some of the few low luminosity sources in 
Fig. E] (z < 0.5 sources from the Sh07 & Mc08 samples) ap- 
pear to lie above our best-fit relation. This might be due to the 
contamination of their optical spectra by host light, which would 
cause their L5100 (and hence Rblr) to be slightly overestimated, 
alth ough effect s hould be very small (< 0.05 dex; see, e.g., Sh07 
and iBentz et al.1 |20o3)). In addition, the few extremely high lu- 
minosity sources (z > 1.5 sources from the D09 & M09 samples) 
also lie slightly above the best-fit relation of the combined dataset. 
Thus, the slope of the .Rblr — ^3000 relation may be somewhat 
shallower or steeper, in the low- or high-luminosity regimes, re- 
spectively. We tested these scenarios by re-fitting the data after 
omitting data points either above L3000 = 10 46 ergs~' or below 
L3000 = 10 45 ergs~', that is keeping most of the (binned) SDSS 
data but omitting the extreme sources from the "small" samples, 
on either side of the luminosity range. The results of this analysis 
are also presented in Table [2] In particular, we find that for sources 
with L3000 > 10 45 ergs _1 the best-fit slope is a ~ 0.7, while for 
sources with L3000 < 10 46 ergs _1 it is a ~ 0.5. Table|2]also lists 
best-fit parameters for other choices of sub-samples. In most of 
these cases, the derived s lopes for the Rblr — ^3 000 relations are 
between those reported bv lMcLure & Jarvisl J2002h (a = 0.47) and 
by MD04 (a = 0.62). The intercepts are also very similar to those 
of iMcLure & Jarvisl J2002h . and the intercept derived from the en- 
tire dataset (1.33) differs from the one derived in MD04 by only 

3 An alternative choice would have been to estimate uncertainties as 
a I yW, Due to the large number of SDSS sources in our sample, this would 
have decreased the uncertainties to below 0.01 dex, which is unrealistic. 
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Figure 4. The relations between Rblr. estimated from L5100, and L3000, for the different "H/3 xMg II" sub-samples under study. Left - large samples: blue 
and magenta points represent the SDSS and 2QZ "HjSxMgll" sub-samples, respectively. Larger symbols and error bars represent the binned data and the 
corresponding scatter (standard deviations). Right - small samples: the Sh07 (red circles), Mc08 (red triangles), D09 (magenta crosses) and M09 (magenta 
triangles) samples. The binned SDSS dataset is also included (blue diamonds and error bars). In both panels, the solid lines represent our best-fit relations 
(Esq. [7}. Other lines provide reference for alternative relations: dashed lines represent our best-fit relations for either high- or low-luminosity sources (slopes of 
0.695 and 0.531, respectively), while the dotted line represents a relation with a slope of 0.5 (normalized to match the median /?blr of the SDSS "Hj8 xMg II" 
sub-sample at L3000 = 10 45 ergs -1 ). See the text and Table[2]for more details. 



0.06 dex, in the sense that our best-fit relation predicts higher ^blr 
values for a given value of L3000 The MD04 study relation was 
derived using only sources with L3000 > 10 44 erg s~ 1 , as is the case 
for most of the sources used here. Therefore, the differences be- 
tween the MD04 relation and our results are not due to different 
luminosity regimes. 

As Table [2] and Fig. [4] demonstrate, the range of slopes and 
intercepts that describe the Rblr — ^3000 relation is relatively 
broad, and probably somewhat luminosity-dependent. This sit- 
uation is present in virtually all previo us attempts to calibrate 
^blr ~ L relations (see discu ssion in, e.g. lKaspi et al .|[200ol . 1 20051 ; 
Vestergaard & Peterson 2006). Notwithstanding this issue, in what 
follows we chose to estimate ^blr by using the relation: 

Ablr feooo) = 21.38 ^ 10 ^g S _! ) It-days. (7) 

We caution that this relation may not be suitable for low-luminosity 
sources (L30OO < 10 45 ergs~'), where a shallower slope should be 
used. 



6.2 Using L (Mgn) to Estimate Rblr 

The line luminosity L(Mgll) can also be used to estimate ^blr- 
Such an approach may be used to overcome the difficulties in de- 
termining L3000 in NIR spectra of high-redshift sourcesQ The rela- 
tion between ^blr> as determined from L5100, and L (Mgll), is pre- 
sented in Figure[5] The scatter in this relation is larger than that of 



4 Several previous studies ha ve calibrated Rblr — L relations for th e Hj3 
and Ha line luminosities (e.g., Greene & Ho 2005; Kaspi et al. 2005). 



the ^blr-^3000 relation. In particular, the 2QZ "HjSxMgll" sub- 
sample shows considerable scatter and almost no correlation be- 
tween Rblr and L(Mgll), probably due to the low-quality of the 
data at this extreme redshift range (see 92.2b . Fitting the data with 
a linear relation of the form 

log^BLR = ! log(^ T A^ T ) +/ 3 (8 ) 

gives the values of a and jS presented in Table |2] We get 
a [L (Mgll)] ~ 0.5 for the SDSS and 2QZ "HjSxMgll" sub- 
samples, but a considerably steeper slope (or a higher intercept) in 
the case we include the high-luminosity sources from the "small" 
samples. The scatter between the resulting (best-fit) estimates of 
^blr an d those based on L5100 (Eq. [2j, for the SDSS and 2QZ 
sub-samples, is about 0.13 dex, only slightly higher than that of the 
^blr-^3000 relations. Despite the advantages of using this Rblr- 
L(Mgll) relation, we draw attention to the significant differences 
in the best-fit parameters that were derived from the different sub- 
samples, as well as by the two fitting methods, and caution that this 
relation is not as robust as the ^blr-^3000 one. In addition, an accu- 
rate determination of L (Mgll) requires a reasonable determination 
of the Fe II and Fe III features adjacent to the Mg II line, and thus 
still depends on the determination of the continuum. 

6.3 The Width of the Mg 11 A 2798 Line 

Our measure of FWHM(Mgll) is different from the ones used 
in several earlier studies that measured the width of the total 
(doublet) profile. At large widths (FWHM(Mgll) > 4000 kins" 1 ), 
the two widths are basically identical, since the width of the 
line is much larger than the separation between the two compo- 
nents. For relatively narrow lines (FWHM (Mgll) < 4000km s" 1 ) 
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Table 2. Coefficients for Rblr-^ Correlations 



Luminosity & sub-samples used 


a 


BCES bisector 
P 


a 1 ' 


M b 


a 


FITEXY 




x = log (L 3 ooo/10 44 ergs" 1 ) 


SDSS 


0.482 ±0.004 


1.440 ±0.004 


0.066 


7.81 


0.531 ±0.004 


1.388 ±0.004 


0.064 


SDSS+2QZ 


0.473 ±0.005 


1.450 ±0.005 


0.077 


7.98 


0.508 ± 0.003 


1.415 ±0.004 


0.076 


"Small" b 


0.619±0.014 


1.340 ±0.022 


0.098 


5.61 


0.620 ±0.0 15 


1.334 ±0.025 


0.098 


SDSS (binned) + "Small" c 


0.618 ±0.014 


1.330 ±0.020 


0.094 


5.62 


0.619±0.014 


1.333 ±0.023 


0.094 


L 30 oo< lO^ergs- 1 


0.531 ±0.015 


1.400 ±0.019 


0.071 


6.93 


0.544 ±0.025 


1.389 ±0.029 


0.071 


/ u;:,,: I0 45 ergs 


0.695 ±0.022 


1.170 ±0.037 


0.087 


4.67 


0.701 ±0.029 


1.158 ±0.059 


0.087 


SDSS+2QZ (binned) + "Small" c 


0.615 ±0.014 


1.340±0.019 


0.093 


5.66 


0.618 ±0.014 


1.330 ±0.023 


0.093 



x = log (i (Mgll) / 1 42 erg s " 1 ) 



SDSS 


0.499 ±0.006 


1.38 ±0.007 


0.124 


0.450 ±0.005 


1.428 ±0.005 


0.123 


SDSS+2QZ 


0.508 ±0.007 


1.36 ±0.008 


0.125 


0.440 ±0.005 


1.436 ±0.005 


0.123 


"Small" d 


0.713 ±0.034 


1.13 ±0.056 


0.151 


0.692 ±0.036 


1.166 ±0.066 


0.151 


SDSS (binned) + "Small" d 


0.684 ±0.032 


1.18±0.051 


0.146 


0.677 ±0.031 


1.190 ±0.057 


0.146 



" Standard deviation of residuals, in dex. 

b The scaling factor associated with the final Mbh estimator, i.e. the equivalents of the factor 5.6 in Eg. 1121 
c The Sh07, Sul07, Mc08, M09 and D09 "H/3 xMg II" sub-samples. 
d The Sh07, Sul07 and M09 "Hj3 xMg II" sub-samples. 
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Figure 5. The relation between i?BLR, estimated from L5100, and L(Mgll), 
for several "H/3xMgII" sub-samples under study. Symbols are identi- 
cal to Fig. [4] The solid black line represent the best-fit relation for the 
SDSS "H/3xMgII" sub-sample, which follows R B LR « L(Mgll) ' 5 . The 
dashed lines alternative relations, derived by including the small, high- 
luminosity and high-redshift sample of M09. These relations follow Rblr 01 
L(Mgll) ' 7 . See the text and Table[2]for more details. 

the two measures differ significantly, and FWHM(MgII, single) < 
FWHM (Mgll, total). For example, the typical profiles with (total) 
FWHM(Mgll) ~ 2000kms~ 1 in our SDSS sample correspond to 
a single-component width of only ~ 1500km s - 1 . This implies that 
Mbh can be systematically overestimated by a factor of ~ 1.8 for 
such narrow-line objects. This issue is crucial for studies of sources 
with high accretion rates. To correct earlier results that used the 



entire profile, we suggest the following simple relation, which is 
based on a fit to the SDSS data: 

FWHM (Mgll, single) = 1.01 x FWHM (Mgll, total) (9) 

-304 kms -1 . 

We used this prescription to correct the relevant tabulated values of 
FWHM(Mg II) published in earlier papers (see TableQ}. 

Several studi es showed that FW HM(Mg II) is very similar to 
FWHM(Hj3) (e.g. IShen et alj|200l and references therein). This 
justifies the use of FWHM(Mgll) as a tracer of virial BLR cloud 
motion. Indeed, the distribution of FWHM (Mgll) /FWHM (Hj8) 
for our SDSS "Hj8 xMg II" sub-sample peaks at -0.05 dex and the 
standard deviation is 0.15 dex. A similar comparison of the IPV 
line widths results in an almost identical distribution. In order to 
further test this issue, we present in Figure [6] a direct comparison 
between FWHM(Mg II) and FWHM(Hj6), for all the "Hj6 xMg II" 
sub-samples. As expected, there is a strong correlation between 
the widths of the two lines. However, in most of the cases with 
FWHM(Hj8) > 4000 km s" 1 , the Mgll line is narrower. This trend 
is also reflected in the binned data shown in Fig. [6] For example, 
for sources with FWHM (HJ3 ) ~ 9000 km s ~ 1 , the typical value for 
the Mg II line is merely FWHM (Mgll) ~ 6000 kms -1 . We find that 
the best-fit relation between these line widths is 

logFWHM(Mgll) = (0.803±0.007)logFWHM(H/3) (10) 

+ (0.628 ±0.027) , 

based on the BCES bisector. This relation is also shown in Fig. [6] 
Our result is in exc ellent agreement, b oth in terms of slope and in- 
tercept, with that o f IWang et al.l ( l2009h . which is based on a subset 
of our SDSS "Hj3xMgIl" sub-sample (about 10% of the sources; 
see 32.11 ). The general trend we find between FWHM(Mg II) and 
FWHM(H/3) is capt ured, in essence, by the st atistical correction 
factors suggested bv lOnken & Kollmeierl I2008T) . which are useful 
for large samples. However, it is not at all clear whether the fit re- 
flects a real, global trend. The alternative is that FWHM (Mgll) ~ 
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Figure 6. A comparison of FWHM(Mg II) and FWHM(Hj3), for the large (left) and small (right) "HjSxMgll" sub-samples. In both panels, symbols are 
identical to Fig. [4] and the dashed lines represent the 1:1 relations. Solid lines represent the best-fit power-law relation, given in Eq. |10| 



FWHM (H/3 ) up to ~ 6000 kms -1 , beyond which there are some 
differences in the mean location of the strongest line emitting gas. 

6.4 Determination of Mbh (Mgll) 

For each source in the "H/3 xMg II" sub-samples we calculated an 
"empirical scaling factor": 

IX (Mgll) = I i? BL R (^3000) FWHM (Mgll) 2 , (11) 

where ^blr(^300o) ls calculated through Eq. [7] In Figure [7] we 
show the distribution of the relevant normalization factor, defined 
as /(Mgll) = M BH (H/3) /ju (Mgll), for the SDSS "HjSxMgll" 
sub-sample. The distribution has a clear peak at a median (mean) 
value is /(Mgll) = 1.33 (1.42), and a standard deviation of 0.32 
dex. This is in agreeme nt with the expected value of / (Mgll) = 
1. For comparison, the iMcLure & Dunlopl ( |2004l) Mbh estima- 
tor was derived assuming /(Mgll) = 0.890 We further investi- 
gated whether /(Mgll) depends on other observables, such as 
source luminosity, SED shape or EW(Mgll). However, we find 
no significant correlations of this type. The only exception, a 
marginal anti-correlation with L/L^dd ( see j )4-2| > is most probably 
driven by the strong dependence of both L/L^dd an d / (Mgll) on 
FWHM(Hj6). In addition, the small scatter in the ^?blr — ^3000 rela- 
tion, in comparison with that of the FWHM (H/3) - FWHM (Mgll) 
relation, suggests that the range of derived /(Mgll) values is 
driven solely by the scatter in FWHM (H/3) /FWHM(Mgll). This 
is supported by the fact that we find no correlation between 
/(Mgll) and ^3000/^5100 while the significant correlation with 
FWHM(Hj6) /FWHM (Mgll) (not shown here) tightly follows a 
power-law with the expected slope of about 2. 



5 The MD04 derivation assumes / = 1 (for Mbh [H/3]), but also introduces 
an offset of about -0.05 dex to minimize the differences between the re- 
sulting MBH(Mgll) estimates and the reverberation results. Our analysis 
considers / (Mgll) to be the factor required to correct for such an offset. 



Next, we combine our best-fit ^blr — ^3000 relation (Eq. [7J 
with FWHM(Mg II) and the median value of / (Mgll) to obtain the 
final form of the Mg II-based Mbh estimator, 



Mbh = 5.6 x 10 6 



£3000 


0.62 


"FWHM(Mgll)" 


2 


lO^ergs- 1 




10 3 kms _1 





This estimator differs from the one given in MD04 in its overall 
normalization, which is higher by a factor of about 1.75 than the 
MD04 one. Thus, we find that the MD04 Mbh estimator causes 
an underestimation of Mbh by ~ . 25 de x. This result is consis- 
tent with the findings of IShen et alj ( l201ll) . where the slope of the 
^blr — £3000 relation was forced to be 0.62. Naturally, the choice 
of a different Rblr — ^3000 relation, according to the parameters 
listed in Table [2] would result in subtle luminosity-dependent dif- 
ferences between our Mbh estimates and those of MD04. We re- 
peated the above steps for several other choices of parameters (a 
and j8 in Eq.[6]l, and list in Table[2]the resulting scaling of the final 
Mbh estimator (i.e. the equivalents of 5.6 in Eq.[l2]above)0 Thus, 
one can use a different L3ooo-based virial Mbh estimator, fully de- 
scribed by a and Mo in Table(2] according to the required luminos- 
ity range (see discussion in 36.11 1. 

To evaluate the improvement we obtained in estimating 
MBn(Mgll), we compare in Figure [8] the masses obtained with 
the H/3 (following Eq. O and the MD04 methods, for all the 
"H/3 xMg II" sub-samples. Clearly, there is a systematic offset be- 
tween the two estimators, in the sense that Mgn(Mg II,MD04) is 
typically lower than M B h(H/3). The median offset within the SDSS 
"HjSxMgll" sub-sample is, indeed, ~ 0.25 dex. In Figure [5] we 
preform a similar comparison, but this time using the new Mg II- 
based Mbh estimates (Eq.ll2t. The systematic shift seen in Fig. [8] 
particularly at the high mass end, has completely disappeared. The 
median difference between the two estimates is negligible, although 
the scatter (standard deviation of residuals) remains about 0.32 dex. 

6 For clarity, we calculated these factors only for the parameters derived 
using the BCES method. 
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Figure 8. A comparison of estimates of Mbh based on the H/J line and the iMcLure & Dunlorl i2004l) prescription, for the large (left) and small (right) 
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Figure 7. The distributions of empirically derived normalization factors 
(/), for different emission lines and sub-samples. The solid red line rep- 
resents the distribution of /(Mgll) (see j|6.4) , for the SDSS "HjSxMgH" 
sub-sample. The solid blue line represents the distribution of /(Civ) for 
the entire SDSS "MgllxClV" sub-sample, while the dashed blue line rep- 
resents the subset of sources which show |v ff(Clv)| < 500 kms~' . 



Several studies suggested to use Mbh estimat ors where the 
exponent of the velocit y term differs from 2 (e.g., I Greene & Hoi 
120051 : I Wang et al .11 20091) . Like any additional degree of freedom, 
it is expected that this approach may reduce the scatter between 
H/3- and Mgll-based estimates of Mbh- The usage of this em- 
pirical approach abandons the fundamental assumption of viri- 
alized BLR dynamics, which is the basis for all mass determi- 
nations considered here. For the sake of completeness, we de- 



rived a relation by combining Eqs [7] and [10] and minimizing 
the systematic offset with respect to the the H/3-based Mbh es_ 
timators. This process resulted in the relation Mbh = 3.62 x 



10' 



£3000 
L 10** ergs- 1 J 



0.62 



j" FW 



111 

2.5 



M . 



The scatter between 



10" ergs-' J [ 10 3 kms- 

this estimator and the one based on H/3, for the SDSS "Hj8 xMg II" 
sub-sample, is of 0.34 dex, almost identical to the one obtained fol- 
lowing Eq. [T2] above. Since this relation departs from the simple 
virial assumption, we chose not to use in what follows. 

Mbh can a l so be obtained by using the Rblr-L (Mgll) relation 
presented in 96.21 We have repeated the steps described above, as- 



suming (R BLR /ld.) ■■ 
and obtained 

M BH = 6.79 x 10 6 



: 24 (L (Mgll) /lO^ergs" 1 ) 



(see Table |2), 



L(Mgn) 


0.5 


"FWHM(Mgll) 


10 42 ergs-!_ 




10 3 kms -1 



1 2 



Mr. 



(13) 

Mbh estimates based on this relation are also consistent with those 
based on H/3, and the scatter for the SDSS "H/3 x Mgll" sub- 
sample is of 0.33 dex, indistinguishable from the one achieved by 
using L3ooo- We note that the Mg II line presents a clear "Baldwin 
effect", i.e. an anti-correlation between EW(Mg II) and L3000 ( e -g- 
lBaldwin|[l977l . BL05). Therefore, the usage of L (Mgll) probably 
incorporates some other, yet-unknown properties of the BLR. 

We tested the consistency of our improved L3ooo-based de- 
terminations of L/L^dd with those based on L5100 (and H/3). For 
this, we calculated L/L^dd f° r the SDSS "H/3xMgIl" sub-sample 
based on the two available approaches: either using the bolomet- 
ric corrections given by Eq.[5]and Mbh estimates given by Eq . 1 121 
or using the bolometric corrections of iMarconi et al ] d2004l) and 
Mbh estimates given by Eq. [3] In all cases we assume Lgdd = 
1.5 x 10 38 (Mbh/Mj;) ergs -1 , appropriate for solar metallicity gas. 
The comparison of the two estimates of L/Lgdd clearly shows that 
the two agree well, with a scatter of about 0.31 dex, and negligible 
systemati c difference. On the oth er hand, using the bolometric cor- 
rection of lRichards et all d2006al 5.15) and the MD04 estimates of 
Mbh results in a systematic overestimation of L/L^dd by a factor of 
about 2.2. 
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Figure 9. A comparison of Hj3-based estimates of Mbh, and the new calibration of a Mgll-based Mbh estimator, for the large (left) and small (right) 
"H/3 xMg II" sub-samples. The H/3- and Mg II-based Mbh estimates are obtained through Eqs.l3land ll2l respectively. In both panels, symbols are identical to 
Fig-E] and the dashed lines represent the the 1:1 relations. 



7 Civ Al 549-BASED ESTIMATES OF M BH 

As discussed in §5\ under the virial assumption, and the known 
ratio of FWHM(C IV)/FWHM(Hj8) in a small number of sources 
with RM-based ^blr measurements, it is possible to use a re- 
liable estimator for the size of the C IV-emitting region (Eq. [4} 
to calibrate a C IV-based estimator for Mbh, by combining 
^blr(^145u), FWHM(C IV) and a known /-factor (/ = 1 in our 
case). Since for the samples where ^blr(C IV) is directly mea- 
sured it is smaller than i?BLR(H/3), by a factor of ~3.7, we expect 
(FWHM (C IV) /FWHM (H/3)) ~ sjfn = 1 .9. In what follows, we 
address the different ingredients of a virial C IV-based estimator, 
and show that for a large number of the type-I AGN studied here 
the C IV measurements are not consistent with the virial assump- 
tion. 

Figures[lO]and[TT]compare FWHM(C IV) with FWHM(Mg II) 
and FWHM(Hj3), respectively. The scatter in both figures is much 
larger than the scatter in the Mg II-H/3 comparison diagram (Fig.[6]l, 
and the line widths do not follow each other. In practical terms, 
for any observed (single) value of FWHM(C IV), the correspond- 
ing values of FWHM(Hj8) covers almost the entire range of ~ 
2,000 — 10,000kms _1 , in contrast with the value expected from 
the assumption of virialized motion and an identical /. We veri- 
fied that alternative measures of the line width, such as the IPV, 
do not reduce the scatter or present any more significant relations 
between the different lines. Most importantly, a significant frac- 
tion of the sources under study (~ 45%) exhibit FWHM (C iv) < 
FWHM (H/3) (Fig. HD, and 26% of the sources in the SDSS, 
2QZ & 2SLAQ "MgllxClV" sub-samples have FWHM(Civ) < 
FWHM (Mgll) (Fig. I lOt. in contrast to the expectations of the virial 
method. The large scatter, and the lack of any correlation between 
FWHM(CIV) and either FWHM(H/3) or F WHM(Mgll), were 
identi fied in several earlier studies of local (e.g. lCorbin & BorosorJ 
Il996t BL05), intermediate- (e.g., S08, F08) and high-redshift (S04, 
N07, Til) samples. We note that the high fraction of sources where 
FWHM(C IV) seems to defy the expectations is not due to a spe- 



cific population of narrow-line sources (i. e., NLSyls), and/or low- 
quality UV spectra, as was suggested by I Vestergaard & Petersonl 
l l2006l) . In particular, we find that 55% of the "H/3 x C IV" sources 
with FWHM(H/3) > 3000kms _1 (i.e., 88 out of 155 broad-line 
sources) present FWHM (Civ) < FWHM (H/3). In addition, 51% 
of the sources in the N07+S04, D09 & M09 samples (33 out of 65 
sources), where the C IV line was measured in high-quality optical 
spectra, also show FWHM (Civ) < FWHM (H/3). 

To further understand the large scatter in Figs. [Tol 
and QT| we tried to look for correlations between either 
FWHM (Civ) /FWHM (Mgll) or FWHM (Civ) /FWHM (H/3) 
and other AGN properties. The different panels in Figure [T2l 
compare FWHM(Civ) /FWHM (Mgll) with L bol , M BH , L/L EM , 
and the shape of the optical-UV SED, av.uvo- I n this com- 
parison, Lbol, Mbh, i/^Edd an d Ofy.uvo were calculated from 
H/3 -related observables whenever possible (i.e. th e M09, Mc08 
and S h07 samples), and using Eq. [3] and the iMarconi et al.l 
d2004l) bolometric corrections. In all other cases these quan- 
tities were calculated from the Mgll-related observables and 
Eqs. [T2] & [5] The shape of the SED is calculated following 
av.uvo = -[log(L 5 ioo/Li45o)/ (5100/1450)] - 1, or a v .uvo = 
- [log (L 300 o/il45o)/ (3000/1450)] - 1, for the "MgllxClV" 
sub-samples that lack L 5 iqo measurements. Most panels show 
considerable scatter, and a lack of any significant correla- 
tions. Although the SDSS sub-sample shows some systematic 
trends of FWHM (C iy) /FW HM (Mgll) with M BH and L/L Edd 
(see also IWang et al.l l201lh . these trends completely disap- 
pear once the smaller samples are taken into account. These 
trends are thus a result of the very limited range in lumi- 
nosity of the SDSS "MgllxClV" sub-sample, and the depen- 
dence of M B h (and L/Ledd) on FWHM(Mgn). In particular, we 
draw attention to the lack of a significant (anti-)correlation be- 
tween FWHM (Civ) /FWHM (Mgll) and a v .uvo- Such an anti- 
correlation is expected due to the similar forms of the ^blr ~ L re- 
lations (Eqs.[2]&[4}, the uniformity of UV-optical SEDs ( jQTj) , and 
the virial assumption. We will return to this point below. The differ- 
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Figure 10. A comparison of FWHM(ClV) and FWHM(Mgll) for the different "MgllxClV" sub-samples. Left - large samples: blue, magenta and red 
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large scatter. 



ent panels in Fig.[l2]also demonstrate that the population of sources 
that do not comply with the basic expectations (i.e. sources with 
FWHM (Civ) < FWHM(Mgll)) cannot be distinguished from the 
general population. 

An alternative way to examine this issue is to drop the as- 



sumption that / is the same for all lines, and consider empirical 
estimates of the /-factors associated with a C IV-based estimator 
of Mgn- T his can be done by combining the Rblr ~ £-1450 relation 
(given by iKaspi et aU^OO! Eq. g) and FWHM(C IV), for all the 
"Hj8 xC IV" and "Mg II xC IV" sub-samples, into a "virial product" 
of the form 



M (Civ) =2.09 x 10 8 



^1450 


0.55 


" FWHM (Civ)" 


2 


lO^ergs- 1 . 




10 3 kms _1 





M . 
(14) 

We can then calculate the C IV-related /-factors (/ (C iv) = 
M^nl ji (Civ)), where Mbh is determined either from L5100 & 
FWHM(HjS) (Eq.0 orL 300 o & FWHM(Mgll) (Eq.EJl. As men- 
tioned above, we expect / (Civ) = /(H/3) = 1. 

The distribution of /(Civ) for the SDSS "MgllxClV" sub- 
sample is shown in Figure [7] It shows a very broad distribution 
with a peak at about /(Civ) = 1.1, in agreement with the ex- 
pected value. Thus, formally, the usage of L1450 and FWHM(C IV) 
can reproduce the correct typ ical Mbh for a large sample of 
sources. The recent study of ICrooml J201 ll) showed that, in- 
deed, for large samples with narrow ranges of redshift and lu- 
minosities, the typical (average) Mbh can b e determined solely 
based on the distribution of luminosities, regardless of FWHM. 
However, the scatter in /(Civ) (the standard deviation), for 
the SDSS "MgllxClV" sub-sample, is about 0.46 dex. More- 
over, 25% of the sources show /(Civ) > 2.1 and an addi- 
tional 25% show /(Civ) < 0.6. This practically prohibits the 
usage of a single scaling factor to determine Mbh in individ- 
ual sources. Since the scatter in £1450/^-3000 an d i-i45o/i-5ioo 
is rather small (see jj4. U , and since the slopes in the different 
^blr ~ ^ relations are similar, most of the scatter in /(Civ) 
is probably due to the scatter in FWHM (C iv) /FWHM (Mgll) 
and FWHM (Civ) /FWHM (H/3). We investigate this point fur- 
ther by plotting, in Figure [T3] the /-factor against several ob- 
served properties. The large scatter, and lack of correlation among 
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samples of different luminosities, is evident in all the panels, and 
in particular in the comparison with a v .uvo- In contrast, a com- 
parison of /(Civ) with FWHM (Civ) /FWHM (Mgll) exhibits a 
very prominent correlation, which closely follows a /(Civ) °= 
(FWHM (Civ) /FWHM (Mgll)) 2 trend. This clearly demonstrates 
that the main source for differences between the different virial 
Mge estimators are the large and unsystematic differences between 
FWHM(C IV) and FWHM(Mg II) (or other velocity measures). 

The recent study by lAssef et aT. I J201ll) presented detailed 
measurements of H/3 and C IV for 12 lensed QSOs at z = 1 .5 — 3.6, 
drawn from the CASTLES survey. These authors suggest that the 
differences between Mbh(H/3) and Mbh(C IV) (estimated using the 



calibration of IVestergaard & Petersonlr2006T) are mainly driven by 
the shape of the UV-optical SED (i.e., L \ 450/^51 00 X and provide 
empirical correction terms to accou nt for these differe nces. As 
noted in 94.11 the (small) sample of lAssef et al. Il l201 J) includes 
heavily reddened sources, with about a half of their sources having 
^5100 > ^1450 (corresponding to a v .uvo < ~ !)• L ess than 15% of 
ou r SDSS sources sh ow such red SEDs. The corrections provided 
bv lAssef et ail ( 1201 l}) thus minimize the effect of the SED shape, or 
reddening, on Mbh estimates. However, as we showed above, the 
main driver for the discrepancy between Mbh(H/3) and Mbh(C IV) 
is most probably related to FWHM(Cl V), and not to the SED 
shape. Indeed, almost half of sources in the lAssef etaT] < l201lh sam- 
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Figure 13. The relations between the scaling factor / (Civ) and several AGN properties (left to right, top to bottom): Lbol, GV.uvo. Mbh and L/L^m. Symbols 
are identical to previous figures. The dashed lines mark /(Civ) = 1, as expected from the virial assumption and /(H/3) = 1 . Note the large scatter in all 
panels, as well as the lack of significant correlations among samples of different luminosities. Here too, the apparent trends within the SDSS "MgllxClV" 
sub-sample (in the lower panels) are due to its limited luminosity range and the dependence of Mbh and L/L^dd on FWHM(Mg II). 



pie have FWHM(Civ) < FWHM(H/3), in contrast with the virial 
assumption, and in agreement with the f raction of such sou rces in 
our "Hj3 xC IV" sub-samples. Moreover. I Assef et alj d201lh report 
a large scatter between FWHM(C IV) and FWHM(Hj8 ), similarly to 
our findings. We also note that the usage of correction terms which 
are based on L5100 (or FWHM(H/3)) is impractical for large sam- 
ples of z ^ 0.8 sources. On the other hand, if the Hj8 line for such 
sources is available (through NIR spectroscopy), it can be used to 
determine Mbh directly (e.g., by using Eq.fSJ, and we see no point 
in using it solely to correct systematics in Mbh(Civ). 

Fi gs . 1 1 2 1 and 1131 demonstrate that the discrepancies associated 
with a C IV-based virial product cannot be accounted for even if 
other basic AGN properties are known. Moreover, some of these 
properties cannot be reliably determined a-priori, given a single 



spectrum that contains only the C IV lineQ We therefore tested the 
possibility of correcting the discrepancy described above by using 
only C IV-related observables, namely the shift (v ff [C IV]) of the 
line center and the equivalent width of the line (EW[ClV]). Pre- 
vious studies suggested that the large blue-shifts often observed in 
C IV indicate that t he emission originate s from non-virialized gas 
motion (e.g., BL05. [Richards et all 201 ll) . Such a scenario may ex- 
plain why a simple virial product (Eq.ll4t fails to scale with Mbh- 
A trend with EW(C IV) might be expected if, for example, there was 



7 While Lbol can be determined from L1450, and a v uvo can be estimated 
given a broad enough spectral coverage, Mbh and thus L/Z-Edd depend on 
the ability of to construct a C IV-based estimator for Mbh- 
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lar offset velocities, which are probably produced by non-virialized 
gas. Thus, C IV cannot be used to precisely estimate Mbh even for 
these specific sources. We also note the great difficulty in robustly 
determining v ff(C IV) in optical spectra of high redshift sources, 
where the systemic redshift determination heavily relies on a few 
UV lines, including C IV itself and the complex, often partially ab- 
sorbed Lya spectral region. 

We conclude that the C IV line is an unreliable probe of the 
kinematics of the BLR gas. There is a significant population of 
type-I AGN, indistinguishable from the general population, for 
which the width of the C IV line contradicts the basic virial ex- 
pectation. A single-epoch Mbh estimator which relies on the width 
of the C IV line provides results that deviate by ±0.46 dex from the 
more reliable H/3-based estimator. This scatter is due to / (Civ) 
only and was estimated assuming a constant /(Hj3) and negligible 
uncertainties in the ^blr — L relations and in the measurements of 
all relevant observables. The uncertainty in indi vidual C IV-based 
estimates of Mbh can be considerably larger (see lWoo et alj < f2oToh 
for the case of/[H/3]). 



Figure 14. The relation between the difference in the Mg II and C IV line 
width, and the offset of the C IV line with respect to the Mg II line, for 
the SDSS "MgllxClV" sub-sample. Larger symbols and error bars rep- 
resent the binned data and the corresponding scatter (standard deviations). 
The dashed line represents FWHM(Clv) = FWHM(Mgll), while the dot- 
dashed line represents FWHM(Clv) = v^FWHM (Mgll), as expected 
from the virial assumption. 



a common origin for the difference in FWHM and the well-known 
Baldwin Effect. 

We first verified that there is no clear trend of nei- 
ther FWHM (Civ) /FWHM (Mgll) nor /(Civ) with EW(ClV). 
In particular, the population of sources with FWHM(Civ) < 
FWHM (Mgll) cannot be distinguished from the rest of the sources 
based on their EW(C IV). Next, We calculated the shifts of the C IV 
line relative to the Mg II for the "Mg IIxClV" sub-samples under 
study. We assumed that the rest-frame center of the Mg II doublet is 
at 2799.1 1 A, and that the C IV profile is centered at 1549.48A. Our 
fitting procedures calculate the observed line centers as the peaks 
of the entire best-fit BLR profile. Figure [14] presents the resulting 
v off (ClV) against FWHM(Civ) /FWHM (Mgll). As is the case in 
Fig. [12] the scatter is very large and covers more than a factor of 
5 in FWHM (Civ) /FWHM (Mgll). However, there appears to be 
a clear trend of an increasing FWHM(Civ) /FWHM (Mgll) with 
increasing C IV f)/«e-shift (negative velocities). 

Figure [14] draws attention to two particular types of sources 
which may, in principle, be used to overcome the problems asso- 
ciated with FWHM(ClV). First, Civ lines in sources with small 
offsets may be dominated by virialized gas. To test this, we se- 
lected a subset of about 4500 sources from the SDSS "Mg II x C I V" 
sub-sample that have |v ff(Clv)| < 500kms~'. The distribution 
of / (C iv) for this subset is shown in Fig. [7] Clearly, the large 
scatter in FWHM (C I v) /FWHM (Mgll), of more than 0.2 dex, 
propagated to the virial products of these sources. Moreover, the 
typical / (Civ) for this subset is 1.49 - considerably higher than 
the expected value (of unity). Second, sources with blueshifts of 
about 1500 — 2500kms~ 1 appear to match the expected value 
of FWHM (Civ) /FWHM (Mgll) ~ 1.9. Here too, the scatter in 
FWHM (Civ) /FWHM (Mgll) for the relevant sources is larger 
than 0.2 dex, and the scatter in / (Civ) is 0.5 dex. Moreover, there 
is no physical motivation for focusing on lines with such particu- 



8 DISCUSSION AND CONCLUSIONS 

The reliable estimators of Mbh an d ^/^Edd obtained in the present 
work should enable us to examine the observed distributions of 
these quantities, and their evolution over cosmic epochs. However, 
since both Mbh ar, d ^A-Edd depend on the source luminosity, their 
observed distributions, at any given redshift, would suffer from se- 
lection effects due to the flux limit of the respective survey. The 
treatment of this issu e requires the app lication of dedicated sta- 
tistical methods (e.g., iKellv et alj|2O10l . NT07), which is beyond 
the scope of this paper. We thus only briefly present here the pre- 
liminary results that concern the upper envelope of the Mbh an d 
L/Z-Edd distributions, and defer the full analysis of the observed 
SMBH evolution to a forthcoming paper. 

We apply the Hj3- and Mgll-based prescriptions to all the 
samples listed in Table Q] For sources with reliable H/3 measure- 
ments (mainly SDSS source s at z < 0.75), the bolo metric luminosi- 
ties are estimated using the iMarconi et al.l d2004l) corrections and 
Mbh ls calculated through Eq. [5] For sources with reliable Mg II 
measurements (mainly SDSS sources at 0.75 ,$ z ~ 2) we use the 
bolometric corrections given by Eq.[5] while Mbh is calculated by 
Eq.[l2] In all cases where more than one line is observed, we prefer 
H/3-based measurements on other measurements. 

Figure Q3] describes the evolution of Mbh with the age of the 
Universe, for the SDSS, 2QZ & 2SLAQ samples (at z £ 2), as well 
as several of the z > 2 samples (N07+S04, M09, D09 and Til). For 
completeness, we also inclu de a small sample of z — 6.2 sources, 
taken from the studies of iKurk et al.1 d2007l) and IWillott et aT] 
d2010h . Fig. Q3] suggests that the most massive BHs grow at the 
fastest rates at z!t 4.5, reach their final masses (Mbh ~ 10 10 M Fi ) 
before z ~ 2, and remain mostly inactive thereafter (see discussion 
in Til). These objects are probably the progenitors of the most 
massive relic SMBHs found in the centers of giant elliptical galax - 
ies (M87 and several other BCGs, see, e.g. lMcConnell et al.l201 lb . 
The most luminous (SDSS) AGN at z ~ 1 are considerably less 
massive, and are probably the descendants of the less luminous 
AGN at z ~ 2, which may be fainter that the 2SLAQ flux limit. 
Such a scenario should be tested critically, by comparing the num- 
ber densities of the different populations of AGN, which trace a 
sequence of increasing Mbh with cosmic time. Finally, Fig.ll5lin- 
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Figure 15. The evolution of Mbh along the age of the Universe for several of the samples studied here: the SDSS, 2QZ and 2SLAQ samples at z, < 2 (small 
blue, magenta and red points, respectively); the combined N07+S04 sample at z — 2.4 and ~ 3.3 (blue triangles); the z > 2 sources from the M0 9 (magenta 
x-sign s) an d D09 (magenta cro sses) samples; the Tl 1 sources at z — 4.8 (black circles); and a sample of z — 6.2 sources, taken from the studies of lKurk et alj 
120071) and lWillott eTaU j2010l black crosses). The M BH estimates are based on Hp (for the N07+S04, D09, M09 and z < 0.75 SDSS sources) or Mg II (for 
the Til, z— 6-2 and z. > 0.75 SDSS sources), following Eqs.l3"1or ll2l respectively. Dotted vertical lines mark z = 1 and 2. 



dicates that the shut-down of accretion onto SMBHs, at z ~ 1—2, 
cannot be associated with a certain "maximal" value of Mbh- 

Figure[T6] shows the evolution of L/L^dd with redshift for the 
SDSS sample studied here0 We draw attention to the smooth tran- 
sition seen in Fig.[T6]at z = 0.75. This region marks the transition 
from Mbh and ^/^Edd estimates that are based on Hj8 and L5100, 
to those based on Mg II and Z^ooo- As explained in 36.41 this tran- 
sition wouldn't appe ar as smooth if these quan tities were obtained 
using th e methods o f iMcLure & Dunlod J2OO4T) and lRichards et all 
d2006tJ) . The steep rise of the upper envelope of the LjL^A dis- 
tribution with redshift (see also NT07) flattens at z ~ 1 , so that the 
highest-L/Lgdd sources approach the expected limit of L/L^dd — 1- 
We stress that the upper envelope of the L/L^dd distribution is not 
affected by selection effects related to the flux limits of the vari- 
ous samples. For example, an un-obscured AGN at z = 1.5, pow- 
ered by a SMBH with M B h = 1O 9 M and accreting at the Ed- 
dington limit could have easily been observed within the SDSS, 
since its observed ;-band magnitude would be ~16.4. Such a source 
would appear as a green point, at z = 1.5 and L/i-Edd = 1 in 
Fig. [16] - a region in parameter space which is clearly dominated 
by sources with much lower Mbh- Thus, Fig. [T6] suggests that 
the vast majority of very massive BHs, with Mbh ~ 10 9 M Q , do 



The 2QZ and 2SLAQ samples add little information in to this figure, and 
will be analyzed in a forthcoming paper. 



not accrete close to their Eddington limit even at z ~ 2, which 
is often considered as "the epoch of peak Quasar activity". Such 
SMBHs have prob ably experienced periods of faster mass accu- 
mulation at z > 3 dWillott et al]|201fi |Pe Rosa et alj|201 ll . Til). 
This scenario is supported by th e extreme rareness of AGNs with 
L hn i > 10 47 atz~ 1 - 2 (see , e.g..lCroom et alj2004| ; |Hasinger et all 



I2OO5I : iRichards et ail l2006bl : iHopkins et alj|200l and references 
therein). As mentioned above, the flux limits have a considerable 
influence on the faintest sources, and practically determine the low- 
end of theobsen'ed distributions we present, at all accessible red- 
shifts. Indeed, deeper surveys (e.g., zCOSMOS, VVDS) have re- 
vealed populations of z ~ 1 — 2 AGN with L/Lpdd — 0-05 (e.g. 
iGavignaud etai]|2008l ; iTrump et alj|2009l ; iMerloni et all2010h . 

To conclude, in this first of two papers we presented a sys- 
tematic study of the methods we use to measure Mbh and i/^Edd 
in type-I AGN at z ~ 1. This was done through a combination of 
larger samples, spanning a broader range in redshift and luminos- 
ity, in comparison to previous works. We focused on the observ- 
ables associated with the Hj3, MgIlA2798 and C IV 11549 emis- 
sion lines, as these are typically available for single epoch spectra 
of high-redshift sources. Our main findings are: 

(i) We provide a luminosity-dependent prescription for obtain- 
ing L[, i from ^3000- The resulting L^ \ is consistent with the one 
obtained from L5100 and differs by a factor of ~ 1.5 from previous 
estimates of this correction, for luminous sources. 
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Figure 16. The evolution of L/Lem with redshift, for the entire SDSS sample studied here. Red, blue and green points highlight sub-samples of sources, 
with Mbh is in the range of 10 7 - 5-7 - 8 , 10 s 5 ~ 8 ' 8 and 10 9 - 5 ~ 9 ' 8 M Q , respectively. Note the lack of very massive BHs that accrete at high rates (i.e. sources with 
M BH > 3 x 1O 9 M and L/L w > 0.2) at z < 2. 



(ii) The width of the Mg II line was shown to follow closely the 
width of the H/3 line up to FWHM ~ 6000km s -1 , beyond which 
the Mg II line width appears to saturate ( 36.3K 

(iii) We obtained an improved Rblr — ^3000 relation, with 
a best-fit slope of 0.62 , consistent with previous studies (e.g., 
iMcLure & Dunlodl2004l) . We further find that the relation is prob- 
ably luminosity-dependent and provide several alternative slopes 
and scalings. Combining this with point (ii) above, we obtained 
an Mbh estimator that is highly consistent with the H/3-based one 
(Eg. 1121). The scatter between the two is 0.32 dex (j j6.4t . This re- 
lation provides Mbh estimates th at are systematically high er, by 
a factor of 1.75, than those of the lMcLure & Dunlopl (2004) M B h 
estimator. 

(iv) The combination of points (i) and (iii) above produces 
Mgll-based estimates of L/L^dd tnat Me systematically lower, by 
a factor of about 2.2, than those derived by previously published 
prescriptions. 

(v) The width of the C IV line shows no correlation with either 
H/3 or Mg II, and for most sources is substantially different from 
the what is expected from the virial assumption. The large scatter 
in the ratios of line widths (~0.5 dex) does not correlate with any 
other observable AGN property, and thus cannot be corrected for 
(SJ7J. The problematic behavior of FWHM(C IV) dominates the as- 
sociated virial products and practically prohibits any reliable mea- 
surements of Mbh using the C IV line. 

(vi) We do not find a significant population of excessively mas- 



sive BHs (Mbh > W 9 M@) mat accrete close to their Eddington 
limit. Such high mass BHs are observed as slowly accreting sources 
(L/L E dd ~ 0.05 — 0. 1), as early as z ~ 2 (Fig. [16}, and might be ob- 
servable in yet deeper surveys. 

A full analysis of the evolutionary trends mentioned here, and oth- 
ers, will be presented in a forthcoming paper. 

Two new relevant papers were published after the submis- 
sion o f the present paper to the journal. The study of IShen & Liu! 
d2012l) presented a new sample of 60 luminous SDSS AGN at z ~ 
1.5 - 2.2, for which the C IV, Mg II, H/3 and Ha lines were mea- 
sured, using high q uality NIR spectroscopy. The conclusions of the 
Sh en & Liul d2012h study regarding the problematic usage of C IV 
as an Mbh e stimator are in exce llent agreement with our findings. 
In particular, IShen & Liul ( 120120 find a large discrepancy between 
FWHM(CIV) and FWHM(H/3), and almost half of their sources 
show FWHM (H/3 ) > FWHM (C iv) . The study of lHo et alj J2012I) 
presented simultaneous rest-frame UV-to-optical spectra of 7 lumi- 
nous AGN at z ~ 1.5, which include the C IV, Mg II, C III] X 1909 
and Ha lines. All Mg Il-based estimates of Mbh were shown to be 
consistent with those based on Ha, while C IV-based estimates of 
Mbh showed large discrepancies, of up to a factor of ~5 ( compared 
with M BH [Ha]). Thus, both the lShen & Liul j2012h and lHo etal.1 
J2012I) studies strengthen our conclusions, using additional high 
quality, high redshift samples. 
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APPENDIX A: PHOTOMETRY-BASED FLUX 
CALIBRATION FOR 2QZ AND 2SLAQ SPECTRA 

Spectra obtained as part of the 2QZ and 2SLAQ surveys, using 
the AAT, are not flux calibrated. Previous studies used single-band 
magnitudes, in either the bj- or g bands, to derive monochromatic 
luminosities, assuming al l spectra follow a uniform UV-to-optical 
SEP of fy oc y- °- 5 (e.g JCroom et al.ll2004l : iRichards et al.ll2005l: 
iFine et al]l2008f ). Here we employed a more robust flux calibra- 
tion scheme, which relies on all the photometric data available in 
the 2QZ and 2SLAQ catalogs. All sources within the 2QZ cata- 
log have tabulated u, bj and r band magnitudes, and all sources 
within the 2SLAQ catalog have the full set of ugriz magnitudes. 
The following two procedures were tested using a unique sample 
of sources, which have duplicate spectroscopy from both the SDSS 
and 2SLAQ surveys, kindly provided by S. Fine. Although the orig- 
inal sample had 639 sources, here we focus only on the 262 sources 
that passed basic quality criteria, and that have z < 1.3, in order to 
avoid the effects of the telluric features on the estimation of Z^ooo- 
The (fully calibrated) SDSS spectra of these sources were fitted us- 
ing our Mg II-fitting procedure (see ij3}. 

First, we test how can the determination of monochromatic lu- 
minosities be improved if one uses the photometric band which is 
as close as possible to the (observed-frame) spectral region under 
question. In the case of L3000, for sources at z — 1, 1.5 or 2, the 
relevant bands are r, i or z, respectively. To test this method for the 
262 sources, we shifted the composite of VdBOl to the appropriate 
redshift and then scaled to the observed SDSS (AB) magnitude in 
the nearest available band. We used the full composite, including 
emission lines, and the scaling was done through synthetic pho- 
tometry, using the entire passband of the relevant filter. Next, we 
measured the continuum flux density of the shifted and scaled tem- 
plate in a narrow band around (rest-frame) 3000A, similarly to the 
main Mg II fitting procedure. The left panel in FigurefATJcompares 
L3000 obtained through this method to the one measured from the 
SDSS spectra. The luminosity estimates are highly consistent with 
those derived from our Mg II fits of the SDSS spectra - , the mean 
difference is 0.06 dex and the scatter around the 1 : 1 relation is 0. 1 1 
dex. 

Since the method described above assumes a uniform SED, 
it necessarily results in calibrated spectra where ratios of XL^ es- 
timated at any two wavelength bands remain constant. In such a 
situation, there is no justification for re-calibrating any specific 
^blr — relation, since all slopes will remain those of the ini- 
tial relation (Rblr~ ^5100 m our case), and the change in scaling 
factors would only reflect the assumed (uniform) SED. Thus, we 
devised an additional method that attempts to overcome this issue. 

We assume that each AAO-observed source can be described 
as a combination of a power-law continuum and an emission (and 
absorption) line spectrum. The power-law component is provided 
by a grid of spectra of the form f v oc v a , with a varying be- 
tween -2 and in steps of 0.05. Each uncalibrated spectrum is 
first smoothed using very broad boxcar averaging (251 pixels, i.e. 
~1080A). Next, the input spectrum is divided by its smoothed ver- 
sion, to establish a (relative) emission line spectrum. The emission 
line spectrum is then multiplied by each of the power-law spectra in 
our grid. All such "combined" spectra are scaled to match the ob- 
served r band magnitude of the source. Finally, the best combined 
spectrum is chosen, such that it best matches the observed g or bj 
magnitudes, for sources from the 2SLAQ or 2QZ catalogs, respec- 
tively. Figure lA2l illustrates the resultant pseudo-calibrated spec- 
trum of J003844.86-004404.5 (z = 0.6074), observed as part of the 
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Figure Al. Comparison of L3000, as measured from SDSS spectra and from our calibrated 2SLAQ spectra, for the 262 sources that have duplicate high-S/N 
spectroscop y of the Mg II spectral reg ion. Left: calibration obtained by the nearest photometric band and assuming the SED follows the uniform composite 
spectrum of lVanden Berk et alj d200ll) . Right: calibration obtained by fitting a power law to all the available SDSS (ugriz) magnitudes. In both panels, blue 
squares indicate luminosities measured directly from the observed spectra, while black crosses represent the results of our Mg II-fitting procedure. The dashed 
lines represent the 1:1 relations. 
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Figure A2. Example of a 2SLAQ spectrum, calibrated using our procedure, 
for object J003844.86-004404.5 (z = 0.6074). The SDSS spectrum is shown 
in black (smoothed by a boxcar filter of 1 1 pixels). The 2SLAQ spectrum, 
calibrated using all the available SDSS photometry and a best-fit power law, 
is shown in red (smoothed over 9 pixels). 

2SLAQ survey, compared with its SDSS spectrum. The right panel 
of Figure [ATI compares Z/jqoo obtained through the power-law fit- 
ting method to the one measured from the SDSS spectra. The mean 
difference between the estimates if now merely 0.03 dex, and the 
scatter around the 1:1 relation is 0.10 dex, indistinguishable from 
the scatter found for the "uniform SED" method outlined above. 
In conclusion, since the two methods provide consistent esti- 



mates of L3000, both can be used in future studies that rely on 2QZ 
and 2SLAQ spectra of type-I AGN. However, only the full power- 
law fitting method allows for a non-uniform UV-optical continuum 
SED, and thus we adopt it in this paper. 



APPENDIX B: ADDITIONAL INFORMATION ON SMALL 
SAMPLES 

The following provides additional information on the literature 
samples used in this work. 

• 59 PG q u asars at 0.25 < z < 0.5, from the study of 
iBaskin & Lao 1 d2005l BL05 in the present paper). For these 
sources, the H/3 line w as observed as part of the survey of 
iBoroson & Greenld 19921) , and fitted using the procedure detailed in 
^3] From the original sa mple of 81 sources, w e omit the 22 sources 
that were also studied in IShang et al] d200l see below). The C IV 
line was observed as part of several HST and WE campaigns, and 
the archival spectra were fitted by BL05. We adopt their C IV mea- 
surements (Table 1 of BL05) but use our own H/3 fits. 

• 44 sources a t z — 2.4 and z — 3.3 from the study of 
dNetzer et al]|2007l N07 in the present paper), which incorporates 
also the sample presented in dShemmer et al. 1 12004| . S04 in the 
present paper). The C IV line was observed in a variety of optical 
surveys and the Hj8 line was observed in a dedicated H- and A"-band 
campaign. All line measurements for this sample were carried out 
using the procedures detailed in ij3] 

• 2 2 PG quasars at 0. 1 < z < 0.4 from the study of dShang et al] 
1 20071 Sh07 in the present paper). The C IV line was observed by 
the HST and the Mg II and H/3 lines were observed in a dedicated, 
almost-simultaneous ground-based optical campaign. This dataset 
is exceptionally useful for our analysis, due to the high quality of 
the observations & related calibrations, and the detailed measure- 
ments of Hj8, Mg II and C IV for all 22 sources. We thus preferred 
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the usage of these measurements over those reported in BL05 (see 
above), for all the Sh07 sources. We used the tabulated parameters 
of the Sh07 study for all lines, as well as the tabulated continuum 
parameters. However, since the Sh07 study fits a power law con- 
tinuum over the entire (rest-frame) UV-optical range, the tabulated 
parameters underestimate L3000 as defined here (see Sj3). We cor- 
rect for this by applying the typical correction factor derived in j(3 

• 90 sources at z < 0.3 from the study of dSulentic et"aLll2007L 
Sul07 in the present paper). This study presents a compilation of 
C IV measurements that are based on archival HST spectra of 1 30 
sources, many of which overlap with the BL05 or Sh07 samples. 
We chose 90 sources for which Hj8 was reliably measured, by 
applying our fitting procedu re to spectra obtained either from the 
SDSS or from the catalog of lMarziani et alj d2003l) . We omitted all 
sources with z < 0.05, to avoid severe host galaxy contamination. 
For the sake of consistency, we pr efer our own fits over the H/3 
line measurements published within lMarziani et al.1 d2003T) . For 11 
sources, the SDSS spectra also cover the Mg II line. These sources 
are thus part of the small subset of sources with measurements of 
all three emission lines. The SDSS counterparts are removed from 

the large SDSS sample. 

• 19 SDSS sources at z ~ 0.36 from the study of dMcGill et al.l 

Eoo| Mc08 in the present paper) for which the Mg II an d H/3 lines 
were observed in a dedicated Keck/LRIS campaign dWoo et all 
120061) . The SDSS spectra of these sources were removed from our 
SDSS sample. 

• 10 sources at z ~ 1 — 2.2 from the study of dDietrich et al.1 
l2009l D09 in the present paper) for which the C IV (9 sources) and 
Mg II (7 sources) lines were observed as part of several optical sur- 
veys. The H/3 line was observed in a dedicated campaign in the J 
and H bands. 

• 3 sources at z ~ 1 — 2.4 from the study of dMarziani et al.1 

M09 in the present paper) for which the Civ line was ob- 
served as part of the Hamburg-ESO survey dWisotzki et alj|200"o)) 
and the H/3 line was observed in the sZ, J and H bands. 

• 4 sources at z — 4.8 from the study of dTrakhtenbrot et all 
EoH Til in the present paper) for which the Civ line was ob- 
served as part of the SDSS and the Mg II lines were observed by a 
dedicated H- and if -band campaign. All line measurements for this 
sample were carried out using the procedures detailed in ^3] 



APPENDIX C: LINE FITTING 
CI The H/3 complex 

We fitted the H/3 line complex in a very similar manner to the 
one described in NT07. The linear continuum is fitted between two 
bands around either 4435A or 4730A and 5 1 IOA. The continuum is 
subtracted from the observed spectrum and an initial fit of the nar- 
row [OlIl]ll4959,5007A lines is used in order to determine the 
systematic redshift and the velocity dispersion of the narrow line 
region (NLR) gas, so to assist with the full H/3-[0 III] fit (see be- 
low). The two [O III] lines are forced to have the same systematic 
shift and width, in the range 300 < FWHM < 1200kms _1 , and a 
flux ratio of 1:3. 

The grid of Fe II template s is constructed by broad ening and 
shifting the original template o flBoroson & Green] Jl992h . The con- 
volution is preformed with a single Gaussian profile with a width 
that increases in steps of 25kms~ 1 , in the range of 1200 to 20000 
km s _1 . The templates are shifted in steps of 69 km s (~lA), in 
the range of -1000 to +1000 kms -1 . Our Fell grid thus consists 



of 23343 different, normalized, templates. The best-fit template is 
chosen through a standard %- minimization scheme, which focuses 
on the flux density in the 4400-4650A range and provides the best- 
fitting scaling factor. The scaled best-fit template is then subtracted 
from the spectra so that the continuum can be re-measured. The 
process is repeated until we receive continuum- and iron-free spec- 
tra. Several of the samples considered here (N07; Sh07; M09; D09) 
were originally fitted by codes that neglected the possible shift of 
the Fell template. However, the differences in the mai n derived 
parameters (L5100 and FWHM(H /3)) are negligible (see iHu et al.1 
l2008l : lMarziani & Sulentidl2012l) . 

Finally, we fit an emission line model which consists of 5 
Gaussians: two broad and one narrow components for the H/3 line 
and a single narrow component for each of the [O III] lines. Each 
of the broad components is allowed to have a width in the range 
1200 < FWHM < 20000 km s" 1 , and to be shifted by as much 
as 1000 kms -1 , relative to the redshift determined by the initial 
[O III] fit. All three narrow components are forced to have the sys- 
temic shift and line widths. These are further constrained to differ 
only by as much as 200 kms~' (shift) and 30% (width), from the 
parameters determined in the initial [Olll] fit, but not to exceed 
the allowed range of 300 < FWHM < 1200kms _1 . The flux ra- 
tio of the two [Olll] features is again fixed at 1:3, while the ratio 
between the narrow H/3 component and the [Offl]A5007 line is 
free to take any value in the range of 0.05 — 1 (i.e. the presence of 
any significant [O III] signal dictates the presence of a narrow H/3 
component). 

C2 The Mg 11 A2798 complex 

We fit a linear pseudo-continuum is to the flux around 2655 and 
3020A. This component completely ignores the Balmer continuum 
emission and thus does not represent the "real" (accretion origi- 
nated) underlying continuum emission. The typical difference be- 
tween the "real" F3000 and the one we measure here is, however, 
small. For example, combining the composite spectrum of VdBOl 
and the best-fit UV-optical power law continuum (/ v oc v~ 044 ) im- 
plies that the measured F3000 is larger by a factor of 1.16. As a 
second test, we used the SDSS sample under study. The "real" 
3000A continuum was estimated by fitting a power law between 
narrow bands at 2200 and 4200A. This limits the test to sources 
with 0.75 < z < 1.2. The observed -F3000 was measured by applying 
the full-scale Mg II-fitting procedure. The resulting median factor 
for these 14532 sources is 1.15 with a standard deviation of 0.06 
dex (10%). We use this factor to convert the "real" L3000 measure- 
ments of the Sh07 samples, to match our L3000 measurements. 

We use a Fe II & Fe III template made of the composite pre- 
pared by Vestergaard & Wilkes (2001) and several additional Fell 
lines in the region not covered by the above composite, kindly pro- 
vided by G. Ferland (private co mmunication). The new lines are 
simila r to the ones presente d in ISigut & Pradha^l d2003l Fig. 13) 
and in lBaldwin et alj d2004, Fig. 5). The Fe II flux under the Mg II 
line is dominated by the red extension of the ~2750A emission 
complex. The additional flux tends to flatten in sources with very 
broad Fe II lines, but the overall effect on the Mg II fitting is very 
small. We consider this template to be more reliable than the 
addition of const ant fl ux under the Mg H line, as done by, e.g., 
KurketalJ d2007) and iFine et all d200£). Oth er studies, such as 
Salviander et alj d2007l) and lShen et alj J2OO8I) , u se templates thai : 
are ve ry similar to ours, following the models in lSigut & Pradhan] 
d2003l) . A visual inspection of the spectra suggests that the iron 
lines blueward and redward of Mg II must not share a common 
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scale factor. Instead, our procedure finds the template which best 
fits the spectrum in the range of 2600-2700A, in terms of shift and 
width of the Iron features. The scaling in the 2900-3030A region 
is allowed to vary by up to 15%. We then combine the two tem- 
plates, characterized by the same shift and width, so that they join 
at2810A. 

Each of the two doublet profiles consists of two broad and 
one narrow emission Gaussians, and a single (narrow) absorption 
Gaussian. We have fixed the intensity ratio of the two broad Mg II 
doublet components to 1:1, suitable for optically-thick lines, while 
for the narrow components the ration is 1:2. Since the narrow emis- 
sion component of Mg II is difficult to disentangle, we have exper- 
imented with an alternative procedure in which the width of the 
narrow Mg II line was bound to match the one measured from the 
[Olll] A5007 line, in sources where both are observed (the SDSS 
"Hj3 xMg II" sub-sample). This technique resulted in either a width 
which is similar to the one achieved by the free-fitting version, or 
in a fit of inferior quality. Our preferred procedure is, therefore, to 
freely fit the narrow Mg II doublet components with the additional 
constraints of 300 < FWHM < 1200 krns -1 , and a shift of no more 
than 200 km s , relative to the systemic redshift. The broad com- 
ponents are allowed to have 1200 < FWHM < 1 5000km s _1 , and 
to be shifted by as much as 1000 kms~' . The (narrow) absorp- 
tion feature is set to have a doublet Gaussian profile with and is 
allowed to have 150 < FWHM < 2000km s~ 1 , and a center which 
is blue-shifted by as much as 5000 km s , with respect to the ex- 
pected (laboratory) wavelengths of the Mg II doublet, assuming the 
systemic redshift. As with the other line-fitting parameters, we ver- 
ified that these choices capture the properties of the vast majority 
of sources that show such absorption features, while not affecting 
the quality of the broad Mg II profile fit. 



Thus, the inclusion of a narrow C IV component is not expected to 
significantly affect measurements of the broad component of the 
line, which we use throughout the present work. This also justifies 
our usage of broad-C IV measurements that include (e.g., BL05) or 
don't include (e.g., Sh07) a narrow C IV component. The narrow 
absorption feature is set in the same manner as in the case of the 
Mg II line (seelC2labove). 

We also assigned a single broad Gaussian component to each 
of the He II A 1640, O III] X 1663 and N IV] X 1486 lines, which 
are forced to have a common line width, in the rage of 1200 < 
FWHM < 8200 kms -1 , and a comm on shift, in the ran ge of -1500 
to +750 kms -1 . As mentioned in iFine et al. these fea- 

tures may in practice account for emission which originates either 
from Fell or other, unresolved, lines (e.g., AlllAl671). The de- 
rived properties of these lines are not used to infer the physical 
parameters of the BLR gas. 



C3 The C IV X 1 549 complex 

Our p rocedure closely follows the one presented in IFine et all 
feoiol) and we refer the reader to the very thorough discussion in 
their Appendix A. In particular, we do not fit the weak Fe II features 
around the C IV line. 

The components of our C IV line model are very similar 
to those of the Mgll line model (see IC2I above). Each of the 
two doublet profiles consists of two broad and one narrow emis- 
sion Gaussians, and a single (narrow) absorption Gaussian, with 
the intensity ratios following the same scaling as for Mgll. The 
broad components are allowed to have line widths in the range 
1200 < FWHM < 15000kms _1 , and to have shifts in the range 
of -3000 to +1500 kms -1 , thus taking into account significantly 
blueshited C IV profi les, which are often observed in type-I AGN 
dRichards et al.ll2oTll and references therein). The narrow compo- 
nents are allowed to have 300 < FWHM < 1200km s _I , and to 
be shifted by as much as 200 kms~' . We note that our choice 
to include a narrow emission component does not mean that such 
a component contributed significantly to the line profile. Indeed, 
several studies claimed that the narrow emission co mponent in 
C IV is very weak, or non-existent altogether (see, e.g., Wills et all 
BL05 and the discussion in IVes tergaard & Peterson! l2006h . 
Although our code allowed the peak intensity of the narrow com- 
ponent to reach 0.5 that of the broad components, its actual con- 
tribution to the overall fitted C IV line profile is indeed found to be 
insignificant. For the SDSS "Mg IIxC IV"sub-sample, the median 
of the EW(C IV,NLR)/EW(CIV,BLR) ratio is only~ 0.007, while 
for 99% of the sources this ratio is smaller than about 0.05 and 
for 25% of the sources it is completely negligible (i.e., < 0.001). 
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